UCLA/06/TEP/27 SLAC-PUB-11936 SPhT-T06/070 hep-ph/0607014 



All One-loop Maximally Helicity Violating Gluonic Amplitudes in 

QCD* 

Carola F. Berger 
Stanford Linear Accelerator Center 

t> ■ 

O ' Stanford University 

o : 

<N . Stanford, CA 94309, USA 

: 

1—5 ■ Zvi Bern 

i> : 

CN . Department of Physics and Astronomy, UCLA 

<N ; Los Angeles, CA 90095-1547, USA 

^ , 

q ! Lance J. Dixon 

^ Stanford Linear Accelerator Center 

'• 

Stanford University 
^ Stanford, CA 94309, USA 

q.: 

Darren Forde and David A. Kosower 
> . Service de Physique Theorique\ CEA-Saclay 

^ \ F -91191 Gif-sur-Yvette cedex, France 

(Dated: July 2006) 



* Research supported in part by the US Department of Energy under contracts DE-FG03-91ER40662 and 
DE-AC02-76SF00515 

t Laboratory of the Direction des Sciences de la Matiere of the Commissariat a VEnergie Atomique of 
France 



1 



Abstract 

We use on-shell recursion relations to compute analytically the one-loop corrections to maximally- 
helicity-violating n-gluon amplitudes in QCD. The cut-containing parts have been computed pre- 
viously; our work supplies the remaining rational parts for these amplitudes, which contain two 
gluons of negative helicity and the rest positive, in an arbitrary color ordering. We also present 

formulae specific to the six-gluon cases with helicities ( — I and (— +H !-+)> as well as 

numerical results for six, seven, and eight gluons. Our construction of the n-gluon amplitudes 
illustrates the relatively modest growth in complexity of the on-shell-recursive calculation as the 
number of external legs increases. These amplitudes add to the growing body of one-loop ampli- 
tudes known for all n, which are useful for studies of general properties of amplitudes, including 
their twistor-space structure. 

PACS numbers: ll.15.Bt, ll.25.Db, 11.25.Tq, 11.55.Bq, 12.38.Bx 
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I. INTRODUCTION 



The forthcoming experimental program at CERN's Large Hadron Collider (LHC) will 
place new demands on theoretical calculations. Finding and understanding new physics in 
this environment will require the study of processes with higher multiplicity than at the 
Tevatron. For example, it is important to improve our understanding of missing transverse 
energy in association with leptons and multiple jets, arising from Standard- Model production 
of W, Z + multi-jets. Such event classes form backgrounds to searches for supersymmetry 
and other models of new electroweak physics. In order to reach the precision required by 
searches for and measurements of new physics, these processes need to be computed to next- 
to- leading order (NLO), which entails the computation of one-loop amplitudes. The crucial 
case of W, Z + 4 jet production — a background to supersymmetry searches when the Z 
decays to a pair of neutrinos — involves the computation of amplitudes with seven external 
particles, including the vector boson. These are challenging calculations. State-of-the-art 
Feynman-diagrammatic computations have only recently reached six-point amplitudes [1, 2, 
3, 4, 5]. 

In this paper, we instead use on-shell methods to compute loop amplitudes. On-shell 
methods were first developed at loop level — in the unitarity method — providing the 
first practical method for obtaining complete amplitudes using previously computed on- 
shell amplitudes [6, 7, 8, 9, 10] 1 . With a few exceptions [8, 13], applications of the unitarity 
method were generally restricted to supersymmetric theories or to the (poly) logarithmic 
part of QCD amplitudes. This limitation arose from the greater complexity of calculations 
using .D-dimensional unitarity [14], required for reconstructing complete QCD amplitudes 
including rational terms. In the past two years a number of related techniques have emerged, 
including the application of maximally-helicity-violating (MHV) vertices [15, 16] to loop 
calculations [17, 18], and the use [19] of the holomorphic anomaly. However, these techniques 
also suffer from the same limitations. Recent improvements to the unitarity method [13, 
20] use complex momenta within generalized unitarity [21, 22, 23], allowing, for example, 

1 Unitarity has of course been a fundamental concept in quantum field theory since its inception (see 
e.g. rcf. [11]). In more recent years, it has become a practical and efficient computational method 
for reconstructing dimensionally regularized [12] amplitudes containing massless particles and multiple 
kinematic invariants. 
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a simple and purely algebraic determination of all box integral coefficients. In ref. [24], 
Britto, Buchbinder, Cachazo and Feng developed efficient techniques for evaluating generic 
one-loop unitarity cuts, by using spinor variables and performing the cut integration via 
residue extraction. Applying these ideas, Britto, Feng and Mastrolia [25] computed the 
cut-containing terms for the most complex six-gluon helicity amplitudes, the terms with a 

scalar circulating in the loop and gluon helicity assignments ( — I 1 h) and ( 1 h+). 

The cut-containing terms for other helicity configurations, and for other components of the 
amplitudes (within a supersymmetric decomposition [26]), were obtained in refs. [6, 7, 24, 
27, 28, 29]. 

The first use of on-shell methods to obtain state-of-the art QCD amplitudes dates to the 
construction of the Z — > 4 parton one- loop matrix elements [21] (or equivalently, by crossing, 
the virtual corrections for pp — > W, Z+2 jets). These matrix elements have been incorporated 
into several numerical programs for NLO corrections to a variety of collider processes [30, 
31, 32]. The technique used in ref. [21] was to first obtain the (poly) logarithmic terms in 
the amplitudes via the unitarity method. The rational terms were found by constructing 
functions with the proper on-shell factorization properties, by matching the known behavior 
as two partons become parallel (collinear singularities), or as invariants formed out of three 
or more partons vanish (multi-particle poles). This rational-term reconstruction method 
did not achieve widespread application because it was unclear how to make it systematic. 
In particular, as the number of external legs n increases, it becomes very difficult to build 
properly factorizing functions; the number of candidate terms in an ansatz for the rational 
terms grows rapidly with n. 

One-loop on-shell recursion relations [33, 34] provide a means for overcoming this diffi- 
culty. They allow for a practical and systematic construction of the rational terms of loop 
amplitudes. Such recursion relations were written down at tree level by Britto, Cachazo and 
Feng [35]. Their work was stimulated by the compact forms of seven- and higher-point tree 
amplitudes [23, 36, 37] that emerged from studying infrared consistency equations [38] for 
one-loop amplitudes (computed using the unitarity-based method), and by the connection 
between twistor space and complexified momenta noted in ref. [39]. Britto, Cachazo, Feng 
and Witten [40] then proved the tree-level on-shell recursion relations using a special con- 
tinuation, or 'shift', of the external momenta by a complex variable z. The proof essentially 
involves only Cauchy's residue theorem, where the residues in the z plane are determined by 
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the amplitudes' factorization properties (for complex momenta). The remarkable generality 
and simplicity of the proof enabled widespread application at tree level [41, 42, 43], including 
to theories with massive particles [44, 45, 46]. 

The extension of on-shell recursive methods to one-loop amplitudes was developed in 
refs. [33, 34, 47, 48, 49]. The essential inputs to this approach are the cut-constructible 
parts of a desired amplitude — those terms containing polylogarithms, logarithms, and 
associated n 2 terms — and appropriate tree and loop amplitudes with fewer legs, which 
appear in factorization limits of the amplitude. Closed-form expressions were found for 
several sequences of one-loop amplitudes with an arbitrary number of gluons. In particular, 
the amplitudes for 'MHV helicity configurations, in which two of the gluons have negative 
helicity and the rest have positive helicity, were computed for the special case that the two 
negative-helicity gluons are adjacent in the color ordering [33, 49]. 

In this paper, we use the same methods to compute the rational parts, and thereby the 
complete expressions, for all remaining n-gluon MHV amplitudes, those for which the two 
negative-helicity gluons are not color-adjacent. Under a supersymmetric decomposition [26], 
n-gluon amplitudes may be thought of as composed of M = 4 and Af = 1 supersymmetric 
pieces together with a non-supersymmetric (Af = 0) scalar loop contribution. Using the 
unitarity method, the M = 4 MHV contributions were computed in ref. [6], and the Af = 1 
terms in ref. [7]. The polylogarithmic contributions to the scalar- loop components of the 
MHV amplitudes were obtained in ref. [18] using MHV vertices [15, 17]; these contributions 
serve as direct input to our recursive construction of the rational terms. 

A number of new aspects must be considered in order to construct loop-level on-shell 
recursion relations. The most obvious feature is that the shifted amplitude A(z) contains 
branch cuts in the z plane. We will handle this feature as in refs. [33, 34], by subtracting the 
cut-containing terms in the amplitude, along with certain additional 'cut completions' that 
cancel spurious singularities, before applying Cauchy's theorem to the remaining rational 
part. 

Secondly, we must ensure that the residues of all poles at finite z are known. The factoriza- 
tion properties of amplitudes with real external momenta are well understood [6, 50, 51, 52]. 
They completely determine the complex factorization properties of tree amplitudes [40]. 
However, the same statement is not true at loop level. In general, non-supersymmetric 
loop amplitudes contain 'unreal poles' which appear for complex momenta but are absent 
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for real momenta. Since these unreal poles are not yet fully understood, the best strategy 
is to consider complex-momentum shifts that avoid channels with unknown factorization 
properties. As discussed in ref. [34], for n-gluon amplitudes the problematic channels always 
have precisely two identical-helicity gluons on one side of the factorization. In this paper, 
we consider MHV amplitudes with two negative-helicity gluons, and we elect to shift the 
momenta of these two gluons. Because these gluons are color-adjacent to positive- helicity 
gluons, the two-particle channels will always have opposite gluon helicities. Hence there are 
no problematic channels in our construction. 

The straightforward application of Cauchy's theorem underlying our computation requires 
that the amplitude fall off as the shift parameter z is taken to infinity. A more elaborate 
construction is needed [34] when this requirement is not satisfied. We can see from known 
four- and five-point one-loop amplitudes that the requirement is indeed satisfied for our 
choices of shift momenta. We assume that the same holds for n > 5. We can confirm this 
assumption at the end of the calculation (and we have done so for n up to 8), by checking 
both factorization properties and reflection or 'flip' symmetries of the computed amplitudes. 

In addition to determining rational terms in amplitudes efficiently, on-shell recursion 
relations can in some cases provide an effective way to compute the (rational) coefficients of 
integral functions appearing in amplitudes [29]. Indeed, this method was used to determine 
the coefficients of all integral functions appearing in the split-helicity n-gluon amplitude, for 
which identical-helicity gluons are nearest neighbors in the color ordering. 

Our interest in constructing all-multiplicity amplitudes stems partly from the desire to 
study the growth in complexity of the amplitudes as the number of external partons in- 
creases. A difficulty with previous computational methods has been the rather rapid growth 
in complexity with the number of external legs. In contrast, with our bootstrap construction 
we find only mild growth, at least for the MHV and split-helicity configurations, allowing us 
to obtain closed-form analytic expressions. Furthermore, experience has shown that analytic 
all-n expressions provide a wealth of intuition into the general structure of the scattering 
amplitudes. 

In a very interesting series of recent papers, Xiao, Yang and Zhu have obtained the 
rational terms for all independent one-loop six-gluon amplitudes using Feynman diagrams, 
by applying spinor simplifications together with integrations that target only the rational 
terms [53, 54]. When combined with the previously-obtained cut parts [7, 18, 24, 25, 27, 36], 



6 



these results provide analytic expressions for all remaining six-gluon helicity configurations 
in QCD. We have compared our results for the MHV case and find complete agreement. We 
have also compared our results for six gluons with the numerical results of Ellis, Giele and 
Zander ighi [4] and find agreement. 

This paper is organized as follows. In the next section we briefly remind the reader of 
our notation. We give an overview of the on-shell bootstrap approach in section III. We 
apply this method to the recursive construction of MHV amplitudes in section IV, giving 
details of the method for five and six external gluons. We also give numerical values at 
selected momentum configurations for the six-, seven- and eight-point amplitudes. These 
values can serve as a reference for future phenomenological implementations. Section V 
contains the all-multiplicity expression for the complete rational parts of one-loop MHV 
amplitudes. We summarize our results and give an outlook in section VI. In the appendices 
we quote previously obtained results for the (poly) logarithmic parts of MHV amplitudes. 
We also give analytical expressions for the newly-computed rational parts of the six-point 
MHV amplitudes. 

II. NOTATION 

In this paper, we shall follow the same notation as defined in section II of our companion 
paper [34]. We therefore refer the reader to that paper for notation. Our amplitudes will 
be expressed in terms of spinor inner products [55]. As is standard for one- loop QCD 
amplitudes, we present color-ordered amplitudes [9, 50, 56, 57]. We need to present only the 
leading-color contributions, because the subleading-color partial amplitudes for a gluon in 
the loop are given by a sum over permutations of the leading-color ones [6]. (For fundamental 
representation particles such as quarks in the loop the leading contributions give the entire 
expression directly.) 

We shall use the supersymmetric decomposition of n-gluon QCD amplitudes [7, 26] to 
rewrite the contributions of different spin states circulating in the loop in terms of super- 
symmetric and non-supersymmetric components, 



A [1/2] _ aM=1 _ 4^=0 
A [l] _ A Af=4 _ aaU=1 



Af=0 
-n;l ? 



(2.1) 



(2.2) 
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where denotes an n-gluon amplitude with a particle of spin J circulating in the loop. 
The non-supersymmetric amplitudes, denoted by M = 0, are just the contributions of a 
complex scalar circulating in the loop, A^° = A^} v The amplitudes A^ 1 and A^ 4 
represent the contribution of an Af = 1 chiral multiplet consisting of a scalar and fermion, 
and an M = 4 multiplet consisting of one gluon, four gluinos and six real scalars, respectively. 

The supersymmetric decomposition is convenient because it separates the amplitudes 
into pieces with differing analytic properties. The supersymmetric parts can be constructed 
completely from unitarity cuts in four dimensions [6, 7], without additional rational contribu- 
tions. The polylogarithms and logarithms of the M = non-supersymmetric contributions 
may also be computed from the four- dimensional unitarity cuts, but in this case there are 
in general additional non-trivial rational contributions. 

The leading-color QCD amplitudes are expressible in terms of the different supersym- 
metric components (2.1)-(2.2) via, 

= <r 4 - 4<r x + (i - es R )A»? + -| (<- - <r°) + ^<r° , (2.3) 

where rif is the number of active quark flavors in QCD. We also allow for a term proportional 
to the number of active fundamental representation scalars n s , which vanishes in QCD. 
We regulate the infrared and ultraviolet divergences of one-loop amplitudes dimensionally. 
The regularization-scheme-dependent parameter 5r specifies the number of helicity states of 
internal gluons to be (2 — e8 R ). For the 't Hooft-Veltman scheme [12] 8 R = 1, while in the 
four-dimensional helicity scheme [7, 58, 59] 8r = 0. 

The amplitudes that we shall present in this paper are unrenormalized. To carry out 
MS renormalization, one should subtract from the leading-color partial amplitudes A n;1 the 
quantity, 



<V 2 e V 3 3 A , 3 .V, 
for D = 4 — 2e dimensions, where, 



n — 2 1/11 2 rif 1 n s 



AT , (2.4) 



i r(i + 6)r 2 (i- e ) 



III. REVIEW OF RECURSIVE BOOTSTRAP APPROACH 



On-shell recursion relations are an efficient way to obtain the rational parts of amplitudes 
directly from their analytic properties. In this section, we first recapitulate the proof [40] of 
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the tree-level recursion relations [35]. The proof relies only on elementary complex analysis 
and on general factorization properties satisfied by any scattering amplitude. The generality 
of the proof permits its extension to loop level, along the lines of refs. [33, 34], as we 
summarize in the rest of the section. 

The proof of the tree-level relations employs a parameter-dependent '[j, /)' shift of two of 
the external massless spinors, j and /, in an n-point process, 

[j, I) : Xj -> Xj - zXi , 

Xi -^Xt + zXj, (3.1) 

where z is a complex parameter. This shift in spinor variables corresponds to shifting two 
of the momenta to complex values, 

a? -fc?w = fc?-|<r|y|r>, 

fcf -^K(z) = k? + z -(r\^\r), (3.2) 

so that they remain massless, k?(z) = kf(z) = 0, and overall momentum conservation is 
maintained. (Similar considerations apply to cases with massive particles [44, 45, 46].) An 
on-shell amplitude containing the momenta kj and k\ then becomes parameter-dependent 
as well, 

A(z) = A(h, . . . , kj(z),k j+1 , . . . , k(z), k l+1 , ...,k n ), (3.3) 

where the physical amplitude is recovered by taking z — 0. 

At tree level, A(z) is an analytic function containing only poles, so we may exploit 
Cauchy's theorem to construct it from its residues. Assuming A(z) — > as z — > oo, then 
there is no 'surface term' in the contour integral around the circle at infinity, so it vanishes, 




Evaluating the integral as a sum of residues allows us to solve for the physical amplitude 
A(0), 

A(0) = - V Res^. (3.5) 

Z Z Oi Z 

poles a 

At tree level, the requirement that A(z) vanishes as z — > oo is rather mild; wide classes of 
shifts exist that satisfy this requirement [35, 40, 44]. 
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FIG. 1: Schematic representation of tree-level recursive contributions. The labels T' refer to 
tree vertices which are on-shell amplitudes. The momenta j and I are shifted, on-shell momenta, 
evaluated according to eqs. (3.14) and (3.15). 

The residues in eq. (3.5) may be obtained using the generic factorization properties that 
any amplitude must satisfy. Consider the factorization channel s Qii ^ — > 0. In this limit, the 
amplitude behaves as [50], 

AT e (h, . . . ,k n ) — > AtTee ( k °" ~ K Lf,) (3-6) 

h=± 

x-xr(v + i,..,ti,ifi), 

S a .../3 

where h — ±1 is the helicity of the intermediate gluon. 

Now consider the effect of the shift (3.2). If the shifted legs j and I are on opposite sides 
of the pole, as depicted in fig. 1, and leg j is between legs a and (5 in the color ordering, 
then under the shift, 

S a ...fi -> Sa...f3 ~ z t a ...n \l~) , (3.7) 

so we have a simple pole in z in this channel. Using elementary complex variable theory we 
may then evaluate the residue, using the rule, 

Res = f{z aP ) . (3.8) 

This rule holds for any analytic f(z) with no additional singularities at z — z a p\ in particular, 
it holds for the product of shifted tree amplitudes appearing as coefficients of each pole in z. 
(Channels with only two isolated legs correspond to a collinear limit, which has a somewhat 
different factorization for real momenta than eq. (3.6). It turns out that for complex 
momenta as needed here, this case can be treated on the same footing as multiparticle 
factorization.) 

Using this rule and summing over all residues then gives us the tree-level on-shell recursion 
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relation, 

AT(h, ...,kn)= J2 E ^ re6 (^i, kp_ lt -K h P ) (3.9) 

partitions P h=± 

x x A tTCC (kp i , . . . , hi, . . . , fcp_ 1 , -^p h ) , 



where 



K P = k Pl + kp 2 + ••• + % + ••• + fcp.! , (3.10) 
K P = k Pl + kp 2 + • • • + kj + ■ ■ ■ + k P _ x , (3.11) 

and h — ±1 labels the helicity of the intermediate state. The hatted momenta are defined 
below. We have defined a partition P to be a set of two or more cyclically-consecutive 
momentum labels containing j, such that the complementary set P consists of two or more 
cyclically-consecutive labels containing /, 

P EE {P 1 ,P 2 ,...,j,...,P_ 1 }, (3.12) 
P EE (P u P 2 ,..., P_!>, 

P U P = {l,2,...,n}, 

which ensures that the sum of momenta in each partition is ^-dependent. The shifted 
complex on-shell momenta kj, ki and K P in eq. (3.9) are determined by solving the on-shell 
condition, Kp = 0, for z, giving the location of the pole, 

" = TOT (3 ' 13) 

Substituting z — > zp in the shifted momenta (3.2) gives, 

k^ ki(zp) = k ^l 1 -J^_ { r\r\r), (3.i4) 

icf = k,(z P ) = k? + \ {n K / pr) (r\Y \r) , (3.15) 

which define kj and k\ appearing in the recursion relation (3.9). 

At loop level, a number of new features arise. In particular, obtaining an on-shell recur- 
sion relation requires dealing with branch cuts, spurious singularities, and in some cases, the 
treatment of factorization using complex momenta, which can differ from 'ordinary' factor- 
ization using real momenta. We refer to channels whose complex factorization properties are 
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non- universal, or at least not yet fully understood, as 'non-standard'. In non-standard fac- 
torization channels, double poles and unreal poles may appear in two-particle channels with 
like-helicity gluons [47, 48] . Fortunately, nonstandard channels for n-gluon amplitudes only 
occur when two adjacent gluons, one of which is shifted and one of which is not, have the 
same helicity. As mentioned in the introduction, for the MHV amplitudes we consider here, 
our electing to shift the momenta of the two negative-helicity legs avoids all nonstandard 
channels. 

To set up a loop- level on-shell recursion we decompose the amplitude into 'pure-cut' and 
'rational' pieces, 

A n (z) = c r [c n (z) + R n (z)\ . (3.16) 

The rational parts R n are defined by setting all logarithms, polylogarithms, and associated 
7r 2 terms to zero, 

R n = -A n =-A n . (3.17) 

Cr rat Cr ln,Li,7r 2 -+0 

The pure-cut part C n consists of the remaining terms, all of which must contain logarithms, 
polylogarithms, or n 2 terms. The cut-containing terms have already been computed for the 
MHV case [18]. Our task is to obtain the rational terms (3.17) via on-shell recursion. 

With the decomposition (3.16), the evaluation of a contour integral like (3.4), but for 
Rn(z), is complicated, in general, by the appearance of spurious singularities. Such singu- 
larities are at kinematical points that do not correspond to any physical singularity. They 
arise in the course of integral reductions and often cannot be removed by algebraic means. 
Although the final amplitude cannot contain unphysical spurious singularities, in many cases 
they cancel only between (poly) logarithmic and rational terms. If we were to set up a contour 
integral (3.4) over the rational terms R n {z) instead of A n (z), it would pick up contributions 
from unphysical poles. 

As a concrete example of a spurious singularity, consider the function, 

ln ( r ) (a-M 

where r is a ratio of momentum invariants. This function has a spurious singularity as 
r — > 1. In the full amplitude this function appears in the combination, 

ln ( r ) + 1 ~ r mQ A 

(l-r-y ' (3 - 19) 

which is finite as r — > 1. 
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To avoid the additional complication of spurious singularities, a good approach is to 
instead 'complete' the cut contribution by replacing functions like eq. (3.18) with non- 
singular combinations like eq. (3.19). Although this procedure is not unique, any function 
that is free of unphysical spurious singularities in z is satisfactory. We denote this cut 
completion by C n , 

6 n (z) = C n (z) + CR n (z), (3.20) 

where CR n (z) are the rational functions added in order to cancel the unphysical spurious 
singularities in z. For a given shift we do not need to remove all spurious singularities, but 
only those that depend on z. Since we have added in rational terms to the completed cuts, 
we should subtract them out from the remaining rational terms. Defining 

R n (z) = R n (z) - CR n (z) , (3.21) 

we have, 

A n (z) = c r [d n (z) + R n (z)j . (3.22) 
With the decomposition (3.22), as shown in ref. [34], the physical amplitude is given by 



A n (0) = Inf A n + c r 



d n (0) - Inf d n + R° + O n 



(3.23) 



where Inf A n is the potential contribution to the amplitude from large z, C n (0) is the 
completed-cut contribution, Inf C n is the potential large- ,2 spurious behavior of the com- 
pleted cut, which must be subtracted off, R® are the recursive diagram contributions, and 
the 'overlap' terms O n remove double counting between the recursive diagrams and the 
rational terms that were added to complete the cuts. 

Let us discuss each of the contributions to the amplitude in eq. (3.23) in turn. First 
consider the Inf A n terms. As discussed in the introduction, using the known four- and five- 
gluon MHV amplitudes [26, 58], one can shift the two negative-helicity gluons and observe 
that the resulting A n (z) falls off as z approaches infinity, that is, Inf A n = 0. We assume 
that this property continues to hold for more than five gluons. We check this assumption 
after the fact, by examining various symmetries and limits of the computed amplitudes. 

Now consider the Inf C n subtraction. In general, since C n (z) is required only to remove 
unphysical singularities at finite values of z, it is possible that the chosen C n (z) has non- 
vanishing behavior as z — > oo. This would introduce spurious large z behavior in R n (z), so 
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FIG. 2: Schematic representation of one-loop recursive contributions. The labels T' and 'L' refer 
to tree and loop vertices. The factorization-function contribution (c) does not appear for MHV 
amplitudes. 

we simply subtract off this large z behavior with Inf C n . If we assume that no logarithms of 
z appear in any surviving terms as z — > oo, the procedure for computing Inf C n is to start 
with C n (z) and series expand around w — 1/z — 0, keeping only the w° term. This term is 
Inf C n . (The assumption can easily be checked in any given case.) In many cases, including 
all cases in this paper, Inf C n goes to a constant at large z, and it can be obtained directly 
by taking the limit, 

Inf C n = lim C n (z) . (3.24) 

z^oo 

The recursive diagram contributions on the right-hand side of eq. (3.23) are obtained 
via a recursion relation, following similar argumentation as at tree level. The sum over 
residues corresponds to a sum over factorization channels. At one loop, as depicted in fig. 2, 
there are in general three contributions to factorization in any given channel s a .„p — > , 

• • • -rl-loop 
^tree x x ^-loop + ^1-loop x ^_ ^ ^toee + ^tree x ^ x ^tree _ 

$a...f3 " Sa.../3 ' s a...(3 

In the first two terms, one of the factorized amplitudes is a one-loop amplitude and the other 
is a tree amplitude. In the Af = scalar loop case, the last term corresponds to a one-loop 
correction to the propagator [51]. Accordingly, in addition to the sum over channels, we will 
have a sum over these additional factorization function contributions. Taking the rational 
parts, we thus obtain, with the shift (3.1), the sum of recursive diagrams, 

- Res ^— = R%(ki,...,kn) 



•2 — Za, 

poles a 



partitions P h=± 



x —j x A tme (k Pi , . . . , ki, . . . , k P i , K p h ) 
+ ^4 trec (^Pi? k P l , —Kp) 

14 




FIG. 3: Diagrammatic representation of overlap contributions. Each overlap diagram corresponds 
to a physical channel. 

, kp_ v K p h ) 
K h P ) 

. . . , fcp_ 1 , K p ) j- , 
(3.26) 

where the partition P is given in eq. (3.12). The 'vertices' R are one- loop amplitudes, but 
with all polylogarithms and n 2 set to vanish, as in eq. (3.17). The term containing the factor- 
ization function JF may be found in ref. [34]. It only contributes in multi-particle channels, 
and only if the tree amplitude contains a pole in that channel. The MHV tree amplitudes 
have no multi-particle poles; hence we encounter no factorization-function contributions in 
this paper. 

The result (3.26) follows directly from the general factorization behavior of one-loop am- 
plitudes, plus the separate factorization of pure-cut and rational terms that was established 
in ref. [33]. Just as in the case of the tree-level recursion (3.9) [40], it exhibits the required 
factorization properties in each channel P (dropping the terms with logarithms, polyloga- 
rithms, and 7r 2 ). Although the R functions are not complete amplitudes, they can be thought 
of as vertices from a diagrammatic perspective, and this equation therefore lends itself to 
the same kind of diagrammatic interpretation available for eq. (3.9). 

Finally, consider the overlap terms O n in eq. (3.23). The recursive diagrams already 
encode the complete residues on the physical poles. Therefore, if the CR n rational terms in 
the completed-cuts C n have any poles in the physical channels, we must subtract these out 
in order to remove this double-count. Since we know CR n explicitly, it is straightforward 
to compute the overlap by performing the shift (3.1) and extracting the residues on the 



x x R(k Pi , . . . , ki, . . . 



+ ^4 tree (^Pn ■ ■ ■ ,kj, . . . , kp_ ± , 

X jp^ X 71 \ K P x i ■ ■ 
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physical poles, i.e., 

poles a 

These overlap contributions may be assigned a diagrammatic interpretation, as depicted in 
fig. 3, with each diagram corresponding to a different physical factorization. Although the 
definition of the completed-cut terms C n is not unique, the ambiguity cancels between C n (0) 
and the overlap terms. 

IV. FIVE-, SIX-, SEVEN-, AND EIGHT-POINT AMPLITUDES 

We will now apply this formalism to the computation of M = (scalar loop) amplitudes 
with two negative-helicity gluons. The logarithmic and polylogarithmic terms C n in the 
N = case are known for all n-gluon amplitudes with two negative-helicity gluons. They 
were computed first for the case where the two negative-helicity gluons are color-adjacent [7], 
and more recently for the general case [18]. We quote the results for the cut parts in 
appendix A, including the specific cut completion we employ in this publication, eq. (Al). 
In appendix A we also explain why this cut completion is suitable: after performing the 
complex shift it contains no z-dependent spurious singularities. 

The rational remainders for all-multiplicity MHV amplitudes were computed previously 
for the special case where the two negative-helicity gluons are adjacent [49]. Below we will 
first illustrate our bootstrap approach by recomputing the rational part of the five-point 
amplitude with two non-adjacent negative helicities [26]. This result also serves as one of 
the inputs for the recursive construction of amplitudes with more than five gluons. Then 
we will compute the remaining, previously unknown, rational terms of all n-gluon M = 
amplitudes with two negative-helicity gluons in arbitrary positions. We start in this section 
with the computation of the six-point amplitudes. In the next section we outline the all-n 
calculation. 

Our bootstrap approach relies on previously-computed amplitudes. A list of amplitudes 
that enter our calculation can be found in the appendix of our companion paper [34]. 
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A. Recomputation of the five-point amplitude A^°(l ,2 + ,3 ,4 + ,5^ 



The five-point amplitude A^ (l~, 2+, 3", 4+ 5+) is given by [26] 

<r°(1^2 + ,3-,4 + ,5 + ) = c r (a 5 + ^ 5 ) , 

where 

C 5 = ^4^(1- 2+ 3" 4+5+) 



(4.1) 



+ i 



(12)(23)(34)(41) 2 [24] 2 2 Ls 



2ro fc.i z^ii -----jJ ^ _|_ Li(~ S23 ) + !,,(" 



-S51 ' — S51 - 



-SB1 ■ 



-S51 ' 



(45>(51>(24> 2 
(32)(2 1)(15)(53) 2 [25] 2 2 Lsi 



'51 



-S12 — S51 ' 
-S34 ' -S34 ' 



+ Li( 



-S12 ' 
-S34 ' 



+ Li( 



-ssi 
-S34 ' 



(54)(43)(25) 2 
2(2 3) 2 (41) 3 [2 4] 3 L 2 



'34 



' -S23 ' 
■ -«51 ' 



3 (4 5) (51) (2 4) 



'51 



2 (21) 2 (5 3) 3 [2 5] 3 L 2 (5gf) 
3 (5 4) (4 3) (2 5) S | 4 



+ 



r 1 (1 3) [2 4] [2 5] ( (1 5) [5 2] (2 3) - (3 4) [4 2] (2 1) ) 



'51 



(4 5) 



2 (1 2) 2 (3 4) 2 (4 1) [2 4f 2 (3 2) 2 (1 5)^(5 3) [2 5] c 
+ 3 (4 5) (51) (2 4) 3 (5 4)(4 3)(2 5) 



+ 



2i 



1 (13) 3 ((15)[5 2](2 3)-(3 4)[4 2](21)) L (5^) 
6 (12) (2 3) (3 4) (4 5) (51) s 51 

(13) 4 2 [2 4] 2 [2 5] 2 



(4.2) 



(12)(41) 2 [24f 



9 (1 2) (2 3) (3 4) (4 5) (5 1) 3 [1 2] [2 3] [3 4] (4 5) [5 1] 3 (4 5) (5 1) (2 4) [2 3] [3 4]s 5l 



1 

+ 77- 



(3 2)(5 3) 2 [2 5]'" 



+ 1 (13) 2 [24][25] 



3 (5 4)(4 3)(2 5)[21][15]s 34 6 s 34 (4 5)s 51 
Here the contribution of an M = 1 chiral multiplet is 



(4.3) 



^(l- 2+3- 4+5+) = 

c r Af c (l-,2 + ,3-,4 + ,5 + 
'(13) 2 (41)[24] 2 Ls 



-K-£-V+ 



+ 2 



(4.4) 



2 T o C "^23 ~^34 
IV-S51' -s 51 



) (13) 2 (53)[25] 2 L Sl ( 



-S12 -ssi 'j 
-S34 ' -S34 ' 



(4 5)(51) 



'51 



(3 4) (4 5) 



'34 



+ 



1 (1 3) 3 ((1 5) [5 2] (2 3) - (3 4) [4 2] (2 1)) L 



-S34 ^ i 

-S51 



(12) (2 3) (3 4) (4 5) (51) 



S51 



We have introduced the function, 
Lsi(ri,r 2 ) 



{l-r 1 -r 2 f 



Li 2 (l - n) + Li 2 (l - r 2 ) + lnri lnr 2 



+ (1 -n -r 2 )(L (ri) + L (r 2 ) 



7T 



(4.5) 
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as well as the L r functions, 




ln(r) 

1 — r ' 




L (r) + 1 

1-r ' 
ln(r) - (r- l/r)/2 
(1-r) 3 



(4.6) 



These functions are free of spurious singularities as r — > 1, by design. 

The task of this subsection is to start with the cut completion (4.2), and use the recursive 
bootstrap to derive the remaining rational terms (4.3). This example has been considered 
already in ref . [33] , but the detailed construction was not given there. We recall the example 
here because the general all-multiplicity solution will follow a similar construction. (For n- 
point amplitudes, however, we will use a cut completion which differs from the above one 
at n — 5, but follows somewhat more naturally from the form of the cut parts obtained in 



We use a [1,3) shift, in the notation of eq. (3.1). The first thing we have to ensure is 



should not produce any ^-dependent denominator factors that vanish at spurious (unphys- 
ical) values of z. Notice that C 5 contains factors of (2 4) and (2 5) in denominators. These 
spurious singularities are not a problem for the [1, 3) shift, because they do not acquire any 
dependence on z. (The double poles in (2 4) and (2 5) are actually cancelled by the behavior 
of the Lsi-containing functions multiplying them.) Spurious denominator factors of the form 
(^23 — s 5i), (si2 — s 34)j and (s 34 — S51) are dealt with via the L 2 functions. 

The next step is to examine the large-parameter behavior of C$(z), which can be obtained 
straightforwardly by applying the [1,3) shift to eq. (4.2). It turns out that in this case, 



The large- ,2 contributions cancel between the series expansion of the (poly) logarithmic terms 
and the rational terms of the cut-completion. This feature is not true in general. Specifically, 
the cut-completion chosen for six or more gluons, as given in appendix A, has a large- z 
contribution (see eq. (A17)). 

As mentioned in the introduction, we 'assume' that the full amplitude has vanishing 



ref. [18].) 



that the cut completion C 5 in eq. (4.2) is a satisfactory one, with respect to this shift. It 




(4.7) 
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large- ,2 behavior under the [1, 3) shift, 

Inf Af= Q = . (4.8) 
[i,3> 5,1 

Of course, in the five-gluon case it is not an assumption, because we know the full amplitude 
in advance. 

With a [1,3) shift, following eq. (3.26), four recursive diagrams give non-vanishing con- 
tributions, as illustrated in fig. 4, 

Rf = Df ] + + Df + Df ] . (4.9) 

As mentioned above, the necessary lower-point amplitudes and vertices are listed in the 
appendix of our companion paper [34], except for i? 4 (l~, 2 + , 3", 4 + ), which is given by, 

R 4 (l~, 2+ 3~, 4 + ) = ( 1 + ^W ee (l-, 2 + , 3", 4 + ) - ^ <13) 2 [1 2] [2 3] 
4V ; V3e 9/ 4 V ; (34)(41)[13] 2 V ' 

(The alert reader will note that Inf[i j3 ) i? 4 = i (1~| 3 1 2 ) 2 / (4 | 1 |3~) 2 ^ 0, due to the last 
term in eq. (4.10). However, this behavior is cancelled by that of the logarithmic terms in 
A^f °, so that Inf [li3 ) A%r° = 0.) 

We obtain for the recursive diagrams, 

D? = Mi~,^, 4+, ^""(S-, - AJ, 2t) = (i + I) (51)(2 g 4)(45) . (4.12) 



= fl 4 (I-,2 + ,A- 3 ;,5 + )— ^~(3-,4 + ,-A-+) 

L + *\ < 13 ) 3 < 14 ) + , < 13 ) 3 (4 13) 

3e + 9j (51)(12)(34)(45)(24) + (1 4) (3 4)(2 5) 2 ' 1 ' J 

W = -^,5 +) ^(s- 4 + , Aa.n = 1 (24>[1 gff 4f [15] • (4.M) 

As remarked earlier, there is no factorization-function contribution for the case of MHV 
amplitudes. 

In addition, there are four channels that can potentially contribute overlap terms, 

[23]' (12)' 

(14)' [35] 1 ; 



% 
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FIG. 4: Non-vanishing recursive diagrams for the amplitude °(1 ,2 + ,3 ,4 + ,5 + ), using a 
[1,3) shift, as given in eqs. (4.11)-(4.14). 

The overlap contributions are illustrated in fig. 5. For our choice of cut completion, eq. (4.2), 
they are given by, 

QW = -1 (12) 2 [24]3 , (4.16) 

3[12][3 4] 2 (4 5)(2 5) V ; 

° (b) = "3 + 2 ) (51) (2 3) (2 4) (4 5) ' (417) 

( C) = » A , 2 \ (13) 3 (14) | i (13) 2 [24](14) .4 (1 3) 2 [2 5](1 5) 



3 Ve / (51) (12) (3 4) (4 5) (2 4) 6 s 34 (45)(l 5) (2 4) 3 s 34 (4 5)(1 4) (2 5) 
,2 (13) 2 [4 5] 2 (13) 3 (4-|(3-5)|2-) , (1 2) 2 (1 3) 3 [23](45) 

% . .. . — T" / j \ / j \/-« . 1 / _ \ Z 



3s 34 (24)(25) 3 534(3 4) (4 5) (14) (2 5) s 34 (l 5)(1 4) 2 (2 4) (2 5) 2 

2 , (13) 2 (15) 2 (24)[25] . (12)(13) 2 [25]((12)(45)-(15)(24)) 

s 34 (45)(14) 2 (25) 2 s 34 (14) 2 (24)(25) 2 

, (13) 3 (15)[35] j ( 13 > 3 f418l 

S3 4(14) 2 (25) 2 (15) (2 4) (3 4) (2 5)' 1 " ) 



i <13>[2 4][3 5]<l-|(2-5)|4-> / (1 4) 2 (1 2)[2 4] 3 



6 [3 4] 2 (4 5)(15)[15](2 4) 3 [2 3] [3 4] (4 5) [1 5](1 5) 2 (2 4) 
i (14)[24] 3 [35] 2 i (15)[4 5] 2 ( gl5 + g35 ) 

3 [2 3] 2 [3 4] 2 (4 5) [15] (2 4) 3 [3 4] 2 (4 5) [1 5] (2 4) (2 5) ' l ' ' 

The final result is the sum of all contributions, eqs. (4.7), (4.11)-(4.14), and (4.16)-(4.19), 

£5(1- 2 + , 3", 4 + , 5+) = - Inf C 5 + D 5 W + L» 5 (b) + ^5 {c) + A> (d) + Of ] + Of ] + 0< c) + Of . 

(4.20) 

The result agrees with eq. (4.3), after appropriate simplification. The complete rational 
part i?5(l~, 2 + , 3~, 4 + , 5 + ), which is needed in the recursive construction of higher-point 
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(c) 



A — 
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l + 5 + 



CR< 




(b) 




C7?< 



(d) 



FIG. 5: Channels giving overlap contributions to A^j~ (l ,2 + ,3 ,4 + ,5 + ) with a [1,3) shift, at 
the values of z given in eq. (4.15). 

amplitudes in the next subsection, is found by adding the rational terms of C 5 and R 5 . It 
can be read off from eqs. (4.2) and (4.3), 



^(l", 2+ 3" 4+5+) = l<r°(l-,2+,3-,4 + ,5 

Cr 



ln,Li2,7r 2 =0 



(4.21) 



B. Six-point amplitudes with two negative helicities 

Without loss of generality, we label one of the negative-helicity gluons as gluon 1 and the 
other as gluon m. Other configurations can be obtained by cyclicly permuting the labels of 
the gluons. We introduce a shorthand for the argument list in the n-gluon case, 



<r°(l, m) = <r°(l-, 2 + , . . . , (m - 1) + , m-,(m + 1)+, . . . , n + ) . 



\Af=0 



(4.22) 



The amplitude where the other negative-helicity gluon is adjacent, A^°(l, 2), can be found 
in refs. [33, 49]. At six points, there are two other independent configurations with two 
negative-helicity gluons, A^=°(l, 3) and A^=°(l, 4). We choose a [1, m) shift, in the notation 
of eq. (3.1), with m = 3,4 in the respective cases. 

The cut parts were obtained in ref. [18]. As explained above, we complete these cut- 
containing terms to remove z-dependent spurious singularities that arise when performing 
the [1, m) shift we've chosen. The result is given in eq. (Al). Taking n = 6 in this expression 
gives the completed-cut parts required here, Cq(0) = Ce(l,m). It is straightforward to 
compute \atC§{\,m) by first shifting Cq(Q) by the [l,m) shift, to obtain C 6 (z), and then 
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series expanding around w — 1/z, keeping the w° term. (The result is given by eq. (A17), 
with n — 6.) 

For the case of Ag=°(l,3) with a [1,3) shift, the recursive contribution is given by, 
itf(l,3) = Af e (-K-,i-,2 + )— J R 5 (3-,4 + ,5 + ,6 + ,K 1 + 2 ) 

+ it! 5 (i-, % 3 , 4+, 5+ 6+)— Af e c (3^, 2+) 

+ Af c (i", X 2 - 34 , 5 + , 6+)— i? 4 (3-, 4+ -K+ 4 , 2+) 

S234 

+ i? 4 (i-, K+ A , 5+ 6+) — e (3-, 4 + , -X- 4 , 2+) 

S234 

+ i? 5 (i-, 2+ X 3 - 4 , 5 + , 6+)— AT e (3", 4 + , 

+ AT c (r, 2 + , X 3 " 45 , 6 +)— J^" 4+ 5 + , -k^) 

s 345 

+ i? 4 (i-, 2 + , i: 3 + 45 , 6+)— Af c (3", 4 + , 5 + , -x 3 - 45 ) 

S345 

+ A 3 ree (i-,-X 6 -i,6 + )— J R 5 (3-,4 + ,5 + ,^ 6 + 1 ,2 + ), (4.23) 

where we have omitted terms that vanish. As mentioned in section III, there is no 
factorization-function contribution for MHV amplitudes. Equation (4.23) is illustrated in 
fig. 6. The amplitudes entering eq. (4.23) are listed in the appendix of ref. [34] and in sec- 
tion IV A. Note that the lower-point amplitudes have at most two negative helicities. This 
general feature of the one-loop MHV amplitudes allows us to solve the recursion relation for 
all n, as shown in section V. 

Once a cut completion is chosen, we can apply the shift [1,3) to CR§ to obtain the 
overlap contributions 0§. For A^f°(l,3) there are six possible channels which can give 
poles, illustrated in fig. 7, 

,(a) = _M ,(b) = M 

[23]' "(12)' 
(c) = s 2M (d) = _ (34) 

(1-|(2 + 4)|3") ' (14)' 

" (1-|(4 + 5)|3-)' "[3 6]' 1 j 

Depending on the specific cut completion, some channels may give vanishing overlap contri- 
butions. 

After assembling and simplifying the recursive and overlap diagrams, we obtain the result 
presented in eq. (B2) of appendix B for i? 6 (l,3), for the specific cut completion (Al). We 
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FIG. 6: Non-vanishing recursive diagrams for the amplitude A^f°(l, 3), using a [1,3) shift, as 
given in eq. (4.23). 




FIG. 7: Channels that can give overlap contributions to °(1,3) with a [1,3) shift, at the 
values of z given in eq. (4.24). 

have checked numerically that our result is symmetric under the flip (123456) <-> (321654), 
even though the recursion relation (4.23) does not keep the flip symmetry manifest. Further- 
more, our result displays the correct multi-particle and collinear factorization limits, as well 
as cancellation of all spurious singularities. We give some numerical values in section IV C. 
The computation of A^f (l,4) uses the [1,4) shift. The recursive contribution is given 

by 

i? 6 D (l,4) = Af c (i-,2 + ,-^r 2 )— J R 5 (4-,5 + ,6 + ,K+,3 + ) 

Sl2 

+ Af»(i- X 2 - 34 , 5 + , 6+)— i? 4 (4", -K+^ 2 + , 3 + ) 
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+ R 4 (l-, K+t, 5 + , 6+)— AT C (4", -k~ u , 2 + , 3+) 

S234 

+ R 5 (i~, 2+ X 3 - 4 , 5 + , 6+)— Af c (4", -^ 3 + 4 , 3+) 

•S34 

+ Af c (i", 2 + , i^, 6+)— it! 4 (4-, 5 + , -X 3 + 45 , 3+) 

s 345 

+ i? 4 (i-, 2 + , ir 3 + 45 , 6+)— Af e (4-, s + , -ir 3 - 45 , 3 + ) 

s 345 

+ R 5 (i~, 2+ 3+ X 4 - 5 , 6+)— Af c (4", 5 + , -X+) 

S45 

+ Af c (i- 2 + , 3 + , X 4 - 56 )— i? 4 (4-, 5+ 6+, -K 4 + 56 ) 

S456 

+ it! 4 (i-, 2+ 3+ X 4 + 56 )— e (4-, 5 + , 6 + , -X 4 - 56 ) 

S456 

+ Af e (i-,-X-,6 + )— ^(4" 5+^+2+ 3+) , (4.25) 

s 61 

as illustrated in fig. 8. Note that the amplitude A^f (l,4) has a flip symmetry under the 
simultaneous exchange of gluons 2^6 and 3^5. This flip symmetry remains manifest 
in the presence of the [1,4) shift. (A second flip symmetry, (123456) <-» (432165), is not 
kept manifest.) Therefore we only need to compute six recursive diagrams in eq. (4.25); the 
other four diagrams can be obtained by the flip symmetry. Diagram (a) can be obtained 
from diagram (j) by the flip symmetry, likewise diagram (b) from (h), (c) from (i), and (d) 
from (g). Diagrams (e) and (f) are invariant under the symmetry. 

With a [1,4) shift there are seven channels that can give contributions to the overlap 
6 (1,4) as illustrated in fig. 9, 



z 



z 



[12] 
[2 4] 

,0 = M 
~ (13) 

(e) = _M 

(15) 

M - M 



(b) _ 




•§234 




(1- 


|(2 + 3) |4-)' 


(d) _ 




S345 




(1- 


1(3 + 5) |4-)' 


(f) _ 




s 456 




(11 


(5 + 6) |4-> ' 



- -[46]- (426) 
Again, we can make use of the aforementioned flip symmetry, 2^6 and 3^5, because 
our cut completion, C7? 6 (l,4) in eq. (A13), was chosen to be flip symmetric. (This property 
relies in part on the symmetry and antisymmetry under r — > 1/r of the functions L 2 (r) and 
Li(r) defined in eqs. (4.6) and (A10), respectively.) The flip symmetry relates overlaps 
and 0( g \ and and and 0^ c \ whereas is invariant. 
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FIG. 8: Non-vanishing recursive diagrams for the amplitude °(1,4), using a [1,4) shift, as 
given in eq. (4.25). 




FIG. 9: Channels giving overlap contributions to °(1, 4) with a [1, 4) shift, at the values of z 
given in eq. (4.26). 

The result for the rational part i2g(l,4), complementing the completed-cut terms given 
in eq. (Al), is listed in appendix B. We have checked numerically that the other, non- 
manifest, flip symmetry (123456) <-> (432165) of this amplitude holds, and that the proper 
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factorization limits are obtained in all channels. We give numerical results for the amplitude 
in the next subsection. 

C. Numerical results for various six-, seven- and eight-point amplitudes 

All-n analytical results were previously found for the M = one-loop amplitudes with all 
gluons of positive helicity [60], and all but one gluon of positive helicity [61]. These results 
have been recomputed via a simplified version of the bootstrap approach [47, 48]. (It is 
simpler because these amplitudes contain no cuts.) The full bootstrap, as described in sec- 
tion III, has been used to compute the n-gluon split-helicity next-to-MHV amplitudes [34], 
using cut-containing terms from ref. [29]. As described in more detail in the next section, we 
have now obtained all-n results for all MHV configurations, by combining the cut parts [18] 
with the rational terms from section V of this paper. The case of MHV amplitudes with two 
nearest neighboring negative- helicity legs in the color ordering was worked out in ref. [49]. 

Here we list some numerical values of the MHV amplitudes for a particular phase-space 
point, defined below, for six, seven, and eight gluons. In the six-point case, we compare 
to results [4] that were obtained with a semi-numerical approach [3], and more recently, 
analytically [54]. For the seven- and eight-gluon cases, we have checked numerically that 
our results have the correct flip symmetries (in fact this was checked up to n = 11), as well 
the correct limits in all multi-particle and collinear factorization channels. 

Specifically, for the six-point case we use the same phase-space point as in refs. [4, 34], 
with the six momenta ki chosen as follows, 

61 = — (— 1, + sin#, + cos6>sin0, + cos#cos0), 

2 

62 = 1, — sin#, — cos sin 0, — cos cos 0), 

k 3 = !(i,i,o,o), 

ki = ^(1, cos (3, sin/3, 0), 

/c 5 = ^ ( 1, cos a cos f3, cos a sin f3, sin a), 

h = —h - k 2 - k 3 - k A - k 5 , (4.27) 

where 

9=j, 0=^, a = |, cos/3 = -^. (4.28) 
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Note that the energies of k\ and k 2 are negative and kf = 0. In order to have energies of 
0(1), the authors of ref. [4] chose the scale n — n — 6 [GeV]. As usual \x also denotes the 
scale which is used to carry the dimensionality of the D-dimensional integrals. 
For the seven-point case we choose the same kinematic point as in ref. [34] , 

k\ = — (— 1, + sin 9, + cos 6 sin 0, + cos 6 cos0), 
2 

k 2 = t^(— 1, — sin#, — cos # sin 0, — cos # cos 0), 

k 3 = |(1, 1,0,0), 

k 4 = — (1, cos/3, sin/?, 0), 
8 

&5 = — (l,cos a cos (3 , cos a sin (3 , sin a), 

k 6 = -^(1, cos (3 cos 7, sin (3 cos 7, sin 7), 

k 7 = —ki — k<i — k 3 — ki — k 5 — /c 6 , (4.29) 

where 

7T 7T 7T 27T 37 , . 

9 — —, 0=—, a = — , 7 = — , cosp = , 4.30 

4 r 6 3 3 128 v ; 

and n = n = 7 GeV. 

At eight points, it turns out that the choice in ref. [34] sits on top of a possible spurious 
singularity (where (7 ± |(1 + 2)|5 ± ) = 0), so here we choose a different reference kinematic 
point, 

h = |(-1, 1/V2, -1/2, 1/2), 

h = -1/V2, 1/2,-1/2), 

k 3 = |(1, 1,0,0), 

A; 4 = — (1, cos/3, sin (3, 0) , 
5 

A; 5 = — (1, cos a cos (3, cos a sin (3, since) , 
6 

/c6 = y (1, cos 7 cos (3, cos 7 sin (3, sin 7) , 

A; 7 = — (1, cos 5 cos /3, cos 5 sin (3, sin 5) , 
8 

ks = —ki — k 2 — k 3 — ki — k 5 — k 6 — k 7 , (4.31) 



where 



8 a 7193 3 3 om 

cosa = — , C08/ ^ = ~8258' cos5 = --, C os7 = -, (4.32) 
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TABLE I: Numerical results for the non-vanishing Af = 4 six-, seven-, and eight-point amplitudes 
with two negative-helicity legs in the FDH scheme. The kinematic points are given in eqs. (4.27), 
(4.29) and (4.31). The analytic expressions used for this table are from ref. [6]. 
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7.274637638 + i 6.314297490 


4.960576389 + i 21.13205060 


-++-++++ 


0.06051809107 + i 0.01883973611 


0.08570866540 + % 0.1142513680 


-+++-+++ 


1.592056562 - i 1.878665237 


5.386800498 - i 1.331984344 



TABLE II: Numerical results for the non-vanishing M = 1 chiral contributions to six-, seven-, and 
eight-point MHV amplitudes. The kinematic points are given in eqs. (4.27), (4.29), and (4.31). 
The analytical expressions were obtained from ref. [7]. 



Helicity 


iA 


e° 


++++ 


- 28.11035867 + i 4.643367883 


- 108.9419206 + i 35.02980993 


-+-+++ 


- 9.102176499 + i 1.503529518 


- 35.86914908 + i 5.750500896 


-++-++ 


0.08664490662 + i 0.6577514371 


- 2.226022769 + i 3.230760457 




- 104.5611840 + i 45.34709475 


- 429.6932951 + i 209.0560823 


-+-++++ 


- 33.85706894 + i 14.68345762 


- 162.2581945 + i 46.94912355 


-++-+++ 


0.6394298622 + i 2.199205804 


- 9.326523676 + i 12.89705727 


++++++ 


- 3.088708769 + i 0.2144746165 


- 7.455275518 + i 0.7776396411 


-+-+++++ 


- 0.3600462794 + i 0.02500099345 


- 0.7753751687 - i 0.06138368438 


-++-++++ 


- 0.002103801109 + i 0.001101616598 


- 0.008916187143 - i 0.003055492662 


-+++-+++ 


0.007480422983 + i 0.09196001806 


-0.2405797542 + i 0.3716908859 
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TABLE III: Numerical results for six-, seven-, and eight-point N = amplitudes as described in 
the text, evaluated at the specific phase space points in eqs. (4.27), (4.29), and (4.31). 



Helicity 


1 1 c 


e 


h+++ 


— 9.370119558 + 1 1.547789294 


— 45.80779561 + 1 13.03695870 


1 III 

— H h++ 


— o.Uo4Uooooo + I U.OUll (DOUoy 


1 K OCtKROOOR _L o' 1 7nOQ7Sn/M 
— ID.ZUDOZZZO + ! 1. 1 Uyo 1 8U44 


i i i i 

— 1 — 1 1 — r 


n 098881 r^^^/i _i_ <; n 9i Q9^n/i7an 
u.uzoooiDo004 -f- % u.ziyzou4/yu 


1 /I /1 08^7007 _L i n 1 7RQ9Q/in^/l 

— i.44yoo/yui -+- 1 u.i ( oozy4U04 


h++++ 


- 34.85372799 + i 15.11569825 


- 176.2169235 + i 87.93931019 


— 1 1 — 1 — 1 — h 


11 98^80^ _L 4 A 80/1/18^879 
— 11 .ZoODoyOO I 4.oy44oDo 1 Z 


^8 977 < ?nfifi/1 _l_ V 1 R7/1' : !889fi 
— Oo.Z / ( OUOD4 -f- % lo.O 1 4000ZD 


-++-+++ 


0.2131432874 + i 0.7330686014 


- 6.238940131 - i 0.4283898751 




- 1.029569590 + i 0.07149153884 


- 4.244770988 + i 0.3284878412 




- 0.1200154265 + i 0.008333664483 


-0.5328645055 - i 0.1225617739 


-++-++++ 


- 0.0007012670363 + i 0.0003672055326 


0.007966818918 - i 0.05081471136 


-+++-+++ 


0.002493474328 + i 0.03065333935 


- 0.1173717691 + i 0.1867703328 



and pL — n — S GeV. 

Our results for the Af = MHV amplitudes are presented in table III. The full QCD 
amplitudes can be reconstructed from the M = amplitudes and the M = 4 and Af = 1 
supersymmetric parts via eq. (2.3). The corresponding supersymmetric amplitudes are listed 
in tables I and II, using results of refs. [6, 7]. Note that we have extracted an overall factor 
of i cr from the numerical values presented in the tables. We do not include the coefficients 
of the leading 1/e 2 singularities of the M = 4 amplitudes in the table, as these are easily 
extracted from the values of tree amplitudes, for any helicity configuration, 



A 



JV=4 



nC T .tree 



(4.33) 

1/e 2 e ^ v i 

The numerical values of the tree amplitudes may be read off from the values of the 1/e 
singularities of either the M = 1 or the M = loop amplitudes, 



—A, 



tree 



jv=o 



AJM = 



C r /.tree 

36 An 



(4.34) 



given in tables II and III. Our results in table III at six points match those of ref . [4] , taking 
into account differing phase conventions. The rational terms in our results agree completely 
with formulae provided by the authors of ref. [54], to high numerical accuracy and for a 
large number of phase-space points. 
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V. ONE-LOOP n-GLUON N = AMPLITUDES WITH TWO NEGATIVE HELIC- 
ITIES 



We now turn to the analytic construction of the all-n amplitude. The complete amplitude 
is decomposed according to, 

<r°(l, m) = c r [4(1, m) + R n (l, m)] , (5.1) 

where 1 and m label the negative-helicity legs. Here C n is the appropriate cut completion 
of C n , given in eq. (Al). To calculate R n (l,m) we will, as in the examples of the previous 
section, choose a [l,m) shift. The remaining rational terms are given by eq. (3.23). Using 
the assumption, extrapolated from the four- and five-gluon cases, that Inf A n vanishes for a 
[l,m) shift 2 , we have, 

R n (l,m) = R°(l,m) + O n (l,m) -IniC n (l,m) . (5.2) 

Our first step is to determine the large- ,2 behavior of the cut-completed terms 4(1, m )- The 
result of this exercise is given in eq. (A17) of appendix A. 

Next we turn our attention to the overlap terms, O n . These terms are calculated by 
taking the residues of the poles in z generated by shifting CR n of eq. (A13) by [1, m). In O n 
we do not include the residues of poles coming from the tree amplitude which is an overall 
prefactor of CR n (l,m) in eq. (A13). Such residues would cancel against corresponding 
terms produced by the insertion of CR n into eq. (5.3) below. Hence we shall drop both 
types of terms from the following discussion. We refer to these modified overlap terms as 
O n . The final remaining piece necessary for calculating the complete rational term comes 
from the recursive rational diagrams, (1, m). These terms are computed using an on-shell 
recursion relation for R n generated by our chosen [l,m) shift. The non-zero terms of the 
recursion relation coming from this shift are, 

R n (l,m) 

= i? n _!(i-, . . . , (m - 2)+, K[ m _ 1)m , (m + 1) + , . . . , n + )—!— 



S(m—l)m 

xAr c (-Kt m _ 1)m ,(m-l) + ,m-) 



2 As mentioned earlier, the assumption that Inf A n vanishes in this case can be checked by requiring 
consistency of the final result with all factorization limits. 
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+ -Rn-l(l )■ 

+ A S (1",2+, 

Here we have combined all the pieces containing known all-n amplitudes [48, 61] into A s , 
which is defined by, 

A s (l-,2 + ,...,m-,...,n + ) 

m min(n,ji+n-3) , 

= E E (^^(1- 2^...,0i-l)^^02 + l)^---,n+)^ 

X % ./ -2i • . . . , m", . . . , j 2 + ) 

+ i^+ii (i-, 2+, . . . , (ji - 1)+, AJ... j2 , (j 2 + 1) + , . . . , n + )-^- 

S jl--32 

x 4 r A 1+ 2(-^.., 2 ^'i + , , J 2 + )) • (5-4) 

As mentioned above, the reason R instead of the full rational term R appears on the right- 
hand side of eq. (5.3), is that the additional CR n terms have been cancelled against cer- 
tain overlap contributions. The three-point loop vertices appearing in eq. (5.4) vanish, 

R 3 (i~, 2+, kU = R 3 (i~, kl {n _ iy n+) = 0. 

The recursive rational contributions naturally split themselves into two classes. The first 
class consists of the one-loop amplitudes with one negative-helicity leg (which are completely 
rational) and the tree amplitudes that multiply them. These terms are fully known and are 
contained in A s . The second class of recursive contributions is, however, more difficult. They 
are given by the first two terms of eq. (5.3) and consist of the rational parts of lower-point 
one-loop amplitudes with two negative-helicity legs, along with the tree amplitudes that 
multiply them. Specifically, these unknown amplitudes are exactly those we are trying to 
determine, having the same number of negative-helicity legs, but with fewer positive-helicity 
legs. The key for obtaining an all-n expression is to solve recursively for this class of terms. 

A. Solving the recursion relation 

To solve eq. (5.3) recursively for R n (l,m), we follow the same general plan as given 
in ref. [46] for the on-shell recursion relations for massive tree-level amplitudes [44, 45]. 



[m 



x AT e (-K 



. , m , n + ) + OJl, m) - Inf CJl, m) . 



(5.3) 
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(Other examples of systematic solutions to on-shell recursion relations may be found in 
refs. [29, 42, 49].) 

Our basic tactic is to insert the left-hand side of eq. (5.3) into the right-hand side of 
eq. (5.3) repeatedly. At each insertion we find that our desired amplitude R n (l,m) appears 
on the right-hand side with one fewer positive-helicity leg, and multiplied by one more three- 
point gluon vertex. This 'unwinding' of the amplitude continues until we have reduced the 
right-hand side of R n (eq. (5.3)) down to, in this case, i? 4 and a sum of terms that contain 
only known quantities (A s , overlap terms O and large-z subtraction terms Inf C) multiplied 
by strings of A^ ee vertices. 

The presence of two terms on the right-hand side of eq. (5.3) containing R n -\ means that 
we remove a leg from either side of m at each step of this unwinding. For example, our 
solution contains the following four terms involving A s at the second unwinding step: 



x^ ree 



A s (r 



xA% ee 



AJr 



xA^ ee 



A s (i~,..., 



xA^ ee 



m - 3)+, K~ 



(m-2)(K (m _ 1)r 



(m + !)+,..., n 



Kt „(m-2) + ,k7 1 



\m-2)(k (m _ 1)m ) 



(5.5) 



(m-2)(K (m _ 1) 

m-2) + ,K 



(m-l)m^ 

,(m + 2)+ . 



■5(m— l)m 



,n 



m-2) + ,K; - ,,(m + 2)- 



Sm(m+1) 



(m-l)(X m(m+1) ) 

^3 ree (-^ + l)^^(^ + l) + )> 

i 



-K+ 



k (m+1) )(. +2 )' K ~^V <•» + 2 ) + )7--^ C (-^ (m+1)! --, (m + D+) 



^m(m+l) 

As can be seen, we get a term for each possible order we can extract two legs from around 
the leg m in the amplitude. Therefore, after extracting / legs (i.e., performing / unwinding 
steps) we obtain 2 l terms. Each term corresponds to a particular ordering of the / legs that 
we extract from around m. 

At each step of the unwinding we must choose new shifted momenta. We always choose 
to shift the two negative-helicity legs of R. For example, after the first step we choose 
[l,i^ m ( m+ i)) or [i, k(m-i)m) as the shifted legs, depending upon which of the two terms on 
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X 2 K 

K[2] 



K m / K xs\ K in 



x 3 




FIG. 10: The propagators corresponding to the sequence of K^y 

the right-hand side of eq. (5.3) we insert i?„_i into. Similarly, when we perform a second 
insertion, of -R„_2, we choose the intermediate K momentum leg of the last shift and the 
previously shifted 1 leg. 

At each step we extract one leg from Ri and create one extra three-point gluon vertex. The 
order of extraction can be thought of as a path of extracted legs denoted by the ordered set of 
momenta {xi,X2, . . . ,Xj + \}, starting with x\ = k m . Then x 2 denotes the momentum of the 
first extracted leg after m, x 3 is the momentum of the leg extracted after X2, and so on, until 
we reach the last extracted leg Xj + i (if j unwinding steps are performed). We need to define 
two sequences of momentum sums related to the chain of intermediate momenta appearing 
as arguments of the three-point gluon vertices. We denote these sequences by {K^}, % = 
1, . . . , j + 1 and {K[{]}, % = 1, . . . , j. We also define a sequence {A^} of corresponding shifts 
of hi, for i = 0, . . . , j. All these momenta are on shell, even The momentum is 

defined recursively by adding the external momentum Xi to -fr^-i] and performing a shift in 
order to bring the sum back on shell. The shift also alters the propagator leg of the previous 
unwinding step from to K[i-i], keeping it on shell. (See fig. 10.) The sequence of 

shifted {k[i\} is dictated by momentum conservation. 

The elements of the sequences involved in the first unwinding step are, for the case 
X2 — m — 1, 



[1] _ h m i 

^m-l)m(l-M fcm) 

2(1 \fi m -i\k ri 
„ K ( m -i)J^\r\k-) 



K » _ U* , Mi _ (m-l)m\ \ I \ ml = " 
-"■[2] ~~ Nn-1 ' K m o/i-IjC l7„-\ — (m-l)m j 



i\ 2 

= k 

2(1 \k m -i\k m 

ftr i — K l , 



K 2 (1 ~\^\k~\ 
"Tii — "'i ~r ; ; ; ; ; — J- 

111 1 2<l-|* m _ 1 |*-> 



(5.6) 
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In this notation, the expression obtained at the first step of the unwinding process is, 

iAf" (-K+(m-l) + , K m ) 

§T " ~Mk m , 2+ . . . , (m - 2)+, K m , (m + 1)+, . . . , n + ) . (5.7) 

(m— l)m 

More generally, after j unwinding steps, we have the following product of A^ ce vertices 
multiplied by an A s term, 



(-1) 



12— m 



"j-j- '• l 3 { n [/•■•'/•■ n i- L] ' 



T=2 



(x r + X[ r -1]) ; 



A s (fc-,, 2 + , . . . , (ix - 1)+, K~ +1] , (i 2 + 1)+, . . . , n+) , 

(5.8) 



where the propagators corresponding to the Ify] are depicted in fig. 10. Here j legs have 
been extracted, (to — i\) legs from the left of to and (i 2 — to) legs from the right of to. The 
factor of (— l)* 2_m , with (^ — to) the number of legs that are after m in the external 
ordering, comes from reversing the order of the arguments of the three-point vertices of legs 
extracted in positions above to in the external ordering, i.e., 

AT{-K+ v K-_ lV x+) = {-l)A^{-K+ v x+,K-_^ . (5.9) 

This reordering is needed so that all extracted three-point vertices have the same order as 
given in eq. (5.8). 

The recursive definitions of the sequences {^[i]} and {%]}, which satisfy the 

on-shell requirements, are 

[1] ~~ m i 

H = < + H-n + «fi , 

H = ^+«^i,. 

*fi = > ( 5 - 10 ) 

where 

(x, + ^- 1] ) 2 (%_ 2] | 7 1^_ 1] > _ (^ + ^_ 1] ) 2 (1-| 7 ^|^-_ 1] ) 
^ 2(^1^1^,) 2(1-1^1^) 

We can make the identification on the right-hand side of eq. (5.11) because the choice of 
shifts is such that the unshifted spinors at step i are given by, 

(*Ml = = • • • = (11 . (^1 = (^[ti]l = --- = (m + \. (5.12) 
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The shifted spinors at step i are given by 



/£+,_/£+ I ■ (zi+i + Kn) 2 lK+] / 7 >-i i K ~\ ^i+i + K^) 2 



(5.13) 



From eq. (5.10) we can see that momentum conservation is satisfied, 

i 

H-i] + H] = +< + H-i] = ■ ■ ■ = k " + E < ■ ( 5 - 14 ) 

Using the representation of the one-loop n-gluon amplitudes with one negative helicity 
in ref. [48], one can see that the left-handed spinor Ai for the negative-helicity leg (leg 1) 
never appears. From eq. (5.12), the right-handed spinor Ai that does appear is unshifted, 

Next, using the identity, 



we rewrite the general term (5.8) as 

(-^-"^(l- xt +1 , xt, . . . , xt, m-) 1 

xA s (r, 2+, . . . , (i, - 1) + , K~ +1] , (i 2 + 1)+, . . . , n + ) . (5.16) 

If A s contains a A m spinor (where m is the negative-helicity leg of A s , other than leg 1) then 
we rewrite the general form of this term as 

(— l) 1 - 2 ^Aj +2 (1 , Xj~ +1 , Xj~, . . . , x 2 , m )—— — 

{lK [j+1] } 2 

xA' s (l~, 2 + , . . . , {h - 1) + , K~ +1] , (i 2 + 1)+, . . . , n + )f(j + 1) . (5.17) 

Here A' s is defined to contain no A m spinors; we construct f(j + 1) such that it contains 
all such factors. The argument (J + 1) denotes the dependence of this function on the 
momentum of the last leg that was extracted for this term in the unwinding, Xj + ±. This 
step is necessitated by our desire to combine together, into a single simple form, all terms 
containing the same set of extracted legs. The terms involved have to be independent of the 
order of extraction of these legs, for it to be possible to combine them. 

Looking at eqs. (5.16) and (5.17) we see immediately that Ky] depends upon the order of 
the path of extracted legs. Because of the homogeneity of the spinors in the amplitude, the 
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denominator factor of (lKy-i) 2 must combine with a complementary factor in the quantity 
A s multiplying it, to generate spinor strings of the form (l~|^ty] . . .. This string can be 
rewritten using (\~\$ {ka b)c = (l~\^ a .. M to get (l"^..^ . . . = (1~\$ h _ i2 . . . (because 



the spinor Ai in the shift always contracts with to give (1 1) = 0). Now, K ix ___ i2 is just 
the consecutive sum of momenta between the left-most extracted leg i\ and the right-most 
extracted leg i 2 . So these terms become independent of the order of the extracted legs in the 
unwinding, and henceforth we will simply replace any such appearance of with K ix ,_, i2 . 

However, this is not the only source of a path dependence. The presence of a A m spinor 
in A s also leads to such a dependence, for example a (ym) term in A s becomes (yKy +1 -\) 
during the unwinding. We find that we cannot completely eliminate the path dependence 
coming from a (ym) term. We can only reduce it to that of the last extracted leg. To see 
this, consider again (yKy +1 j) as it will appear in the amplitude combined with (lKy +1 j) , 



homogeneity of the spinors. From this example we see that, after removing any dependence 
upon the shifted external momentum from the amplitude, we can distinguish different paths 
of extracted legs only by the last leg extracted from R. For the amplitudes we consider here 
we find that such X m spinors are always present, so to proceed we must first isolate all X m as 
in eq. (5.17). Once they are isolated, we can straightforwardly reduce the path-dependence 
to that of the final extracted leg only. 

With this simplified dependence upon the extracted path we can now proceed to combine 
together all terms containing the same set of extracted legs {ii, . . . , to, . . . , i 2 }. All terms that 
correspond to each possible path between two legs are combined. For example, considering 
the case when i ± — to — 1 and i 2 = m + 1, we combine the middle two terms of eq. (5.5), as 
these correspond to the two possible ways of extracting the set of legs {m — 1, m, m + 1}. In 
general for the set of extracted legs {ii, . . . , i 2 } we see that the last leg extracted from such 
a contributing term can only be either ii or i 2 . Therefore we write the sum of all possible 
extraction paths of the legs {i±, . . . , to, . . . , i 2 } as two sums, one for Xj + i = k ix and one for 
x j+ i = k i2 , to get, 





5f,<7(a;i,/3i),TO ) 
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<reOP(a 2 ,/3 2 ) 



x — -4(1" 2+ > • • • > (<i " (<a + !) + ' • • • ^ + ) 

< 1 ^ii...i 2 > 2 



CTeOP(a 2 ,/3 2 ) 

2 



X ■ 



2 4(1- 2+ . . . , (t! - 1) + , K-^ i2 , (i 2 + 1)+, . . . , n + ) , (5.19) 



(l«ii...«> s 

where cr(o!j, /%) indicates one of the possible 'ordered permutations (OP)' or 'mergings' of the 
ordered sets a, and Pi, with = {«i+l, . . . , to— 1}, /3i = {i 2 , . . . , m+1}, a 2 = . . . , m— 1} 
and /9 2 = {i 2 — 1, . . . , to + 1}. The elements of each merging are the union of the two sets, 
ctj U Pi = {xj, . . . ,x 2 }. The merged ordering must preserve the order of the elements ai 
and Pi individually, but any relative ordering of elements of a, with respect to those of Pi is 
permitted. Also, np i is the number of elements in the set Pi. We can then rewrite eq. (5.19) 
using the Kleiss-Kuijf relation [62] between tree amplitudes (for which a simple proof was 
given in ref. [63]). The Kleiss-Kuijf relation reads, in our notation, 

AT(l,a,n,p T ) = (-l)^ AT(l,°(a,P),n), (5.20) 

cr60P(a,/3) 

where P T is P written in the reversed order. Applying it to eq. (5.19) and rearranging the 
resulting MHV tree amplitudes to restore gluon 1, we obtain, 



(t2l)(iim)/(ti) - (ii \){i 2 m)f(i 2 ) AtICC , + _ 

(lm)( ?2ll ) A 2 -h + 2(1 •••,"» ,...,i 2 ) 



x ' 4(1~, 2+, . . . , (i, - 1)+, kr ( i2 + 1)+, . . . , „+) . (5.21) 

^4..',) 

The final result is then given by summing eq. (5.21) over all possible sets of extracted legs 
{ii, . . . , to, . . . ,i 2 }, which is equivalent to summing over all possible factorization channels 
of the amplitude. 

In the simplest case, A s contains no \ m spinors and f(i) = 1. In this case, using the 
Schouten identity, eq. (5.21) reduces to 

^i 2 4+2(l ,iiT--i m j---j*2~) 

X 71 * X2 Ml~, 2 + , . . . , (h - 1)+, tfr ( ia + 1)+, . . . , n + ) . (5.22) 
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The case of the amplitude considered here is more complicated. We always have to take 



into account a \ m spinor contained in the denominator factor of s 



j 1 ...m...j 2 



in the A* of 



eq. (5.4). During the unwinding, Sj 1 ...m,...j 2 becomes Sj 1 ,,,K [j+1] ...j 2 y which, using eq. (5.18), can 
be rewritten as 



Sj 1 ...K b+1] ...j 2 — (^■ + i]|-^ji...(u-l),(i 2 +l)...j 2 |-^[j+l]) + S ji--(h-l)(i'2+l)-j2 



(X 



where we define, 



(x j+1 l) 



a...(b— l),(c+l)...d i 
G a ,b,c,d = fta...(b-l),(c+l)...dfta...d- 



(5.23) 



(5.24) 



and again Xj+i is the last leg extracted during the unwinding. 

Hence, if we isolate a factor of l/s_j 1 ...m...j 2 fr° m the remainder of A s (along the lines of 
eq. (5.17)), then we see that f(j + 1) for the final extracted leg Xj +i is given by 

(x j+1 1) 



f(j + 1) 



(5.25) 



( X j+ll^ji,ii,i 2 J 2 |l + ) 

Inserting this result for f(j + 1) into eq. (5.21) and taking into account the cases when 
i\ = m or i 2 = m we get the result for a specific set of extracted legs . . . , m, . . . , i 2 }, 



i x min(n,n-3+ji) 

-E E 



Ji=2 j 2 =max(i 2 ,i 2 +2-ii+ji) 

x A'J^(1-, 2 + , . . . , (*! - 1)+, X- ., 2 , (i 2 + 1)+, . . . V 



1 K Y- 

- 1 - ±l ll...l2/ 



(5.26) 



with 



JC(a,b,i 1 ,i 2 ,m) 



11 = m and i 2 = m : 1/ s a 5 , 

<! = "»: (<2l>/(<2-|Si il>i3 ^|l + > , 

12 = m: (Hi)/ (l~\Ga,h,i 2 ,b\ii + ) , 

otherwise : 



(5.27) 



("»-|^i 1>ia> Jl + ><l<l><<2l> 



(ml)(l-|^ 1 , l2 , fe K 1 + )(^|^ iil2i6 |l + ) 



and 



A^ 2 (l-,2 + ,...,m-,...,n H 
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= 2 + , . . . , Vi ~ K-.., 2 , (h + 1) + , • • • , n + ) 

x /»';, r .2! A'; ... /V ./;. . . . , m~, . . . , it) 

+ R n . n+n (i- 2+ . . . , (ji - 1) + , AJ. .. j2 , (j 2 + 1) + , . . . , n+) 

x 4^ 1+2 (-i^.., 2 , A + , . . . ,m" . . . ,j 2 + ) , (5.28) 

which are the terms of A s with a factor of « and the denominator factors of s- ■ 
removed. 

To obtain the final result we must include all possible factorization channels. We also 
must include the overlap and Inf C n terms. Finally, there is a contribution from the terminal 
step of the unwinding, because R4 is non-zero for this amplitude. The complete result for 
R n is then given by, 

R n (l,m) 

m min(n,n+ii-5) 

= "E E ^ Ii+2 (r,^...,m-...,i^ 

U=2 i2=m 

I Inf 4-^(1- 2+ . . . , (2! - 1) + , K' ^, (i 2 + 1) + , . . . , 71+) 
min(n,n— 3+ji) 

+ E E fc(ji,j2,ii,i2,m) 

ji=2 j 2 =max(i 2 ,j2+2-ii+ji) 

x 2+, . . . , (ij - 1) + , X- ., 2 , (i 2 + 1) + , . . . , n + ) 

+ 0^g +n (l-, 2+ . . . , (<! - 1) + , ., 2 , (i 2 + 1) + , . . . , n+)) 
- ^-(1-, 3 + , . . . , (m - 1) + , m" (m + 1) + , . . . , (n - 1' ! 



(l^u...i 2 > 2 



(n2) 2 

+ |i4*-(l,m). (5.29) 

Here A'/ 1 - 72 is as given in eq. (5.28) and in a corresponding treatment 0' jlj2 represents 
the terms of O containing denominator factors of s - . i .■„, but with those factors 

extracted; as in A'/ 1 - 72 we also extract an overall factor of i. Furthermore, as mentioned 
above, O- 71 - 72 does not include terms that stem from the residues of poles of the overall tree 
term in CR n , because they are precisely cancelled by CR n . For this reason, CR n is not 
included in eq. (5.29). Finally, note that for i± — i 2 — m, the 'tree amplitude' in eq. (5.29) 
is to be evaluated as A 2 rcc (l~, m~) = i(l m) 4 /((l m)(m 1)) = — i(l m) 2 . 
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B. The final result 



After combining all the necessary terms into eq. (5.29) we find that the final result for 
R n is given by, 

R n (l,m) 



x 



x 



m min(n,n+ii-5) , i 1 min(n,n-3+ji) ,. . . . . .,. i \ • \ / • / • , i \\ 

(r V ^J^J^h^rn) {{ii - l)ii){i2 {i2 + 1)) 

2^ 2^ 2^ 3 (Hi)(i 2 l) 

il=2 i 2 =m x ji=2 j 2 =max(i 2 ,i2+2-ji+ii) N ' N ' 

Knn Hggj; Wggj; W + g^>) + 4...(n-i),fe+i).., 2 (^u 2 H + g^>))) 
(O'l - 1)1^,^ |1 + ) (l-|^Ln,. 2 J 2 I0'2 + + ) 



+ T 6 ill2 (m,n) +T 7 ^ 2 (n)J + T 5 (m,n) j , (5.30) 
where we have introduced 

Ga,b,c4 — fia...(b-l),(c+l)...dfib...c, 

Ga,b,c,d = ftb...cfta...(b-l),(c+l)...d, (5.31) 

in addition to the definitions in eq. (5.24). This solution then depends on the functions, 

ji = h '■ l/ s (i 2 +i)...j 2 , 

= { h<il . <0i-1)Ji> <l-|g|l + > (5-32) 
(l-|^b-i + )((ii-l)-|^t|l+)' 



32 = h ■ l/Sji...^-!) , 
^» = \ . ^ . (j 2 (j 2 +l)>(l-|£|l + > ( 5 - 33 ) 



32 > %2 ■ 



2 -|S|l + }<l-|S|(i 2 + l)- 



where above and throughout this section Q = Gj u i u i 2 j 2 , = <3j u i u i 2t j 2 , Q = GhMM,h and 



Q ] = Q\, h j, „• The function T x is given by, 



J1,*1,»2,J2" 

j 2 ii-2 



;=i 2 +i i=ji-i 



(5.34) 
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with 



ji<l<k-2: 

Q$ l{l+1) $ [l+1) ... {n ^ {l-\Q (11 G \{l + 
(1(1 + 1)) 

I = and j 1 =i 1 : , 

/ = ji - 1 and ji < h : 

-Sji...(»i-i),(i2+i)...j2 ^0i+i)...(u-i),(i2+i)...j 2 ^i...(n-i)^ t + ) 

(i-|^,., 2 ^,., 2 |i + >(i-|^b'i + ) 



(l-ls f b'l + ) 



(5.35) 



i = J2 : 

s ji...(ii-l),(i2+l)...j 2 

(l-|^b 2 + )(l-|^ 1 ... j2 ^ 1 .., 2 |l + ) 



0' 2 -|^|i + ) 

i 2 + 1 < Z < j2 ■ 

(i-\G$i(i + i)$( i2+ i)...i& <i-|a n (i-\g\(i + 1)- 



(/(/ + !)) 



Next, the T 2 term reads, 



n 



h-3 



= E /S 2 (^)+ E 

;=j 2 +2 /=ii— i 



(5.36) 
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with 



ji < I < (ii ~ 1) : 



_ ((/-!)/) (11 g|t L „ p ^ (p+1) ... fa _ 1) gt|i + ) 3 

(P(P+ 1)> 01 ^...(n-l)^^)] 2 ^^!)...^-!)^ + ) 

x s ,.. p (i-| \ P +) (i-| Gti... i i 1 -i ) ti... P \ip + 1) + ) ' 

Ji : 

(11^ IJi + ) (11 ^...^^...(n-i)^ Il + ) 3 



_ y .. 

p4l+i I 0#(P+i)-(ii-i)#ii-i£ (1 I Iji + ) 

(p(p+ i)> (ii g^ 1 ...fa-i)TOi,p)] 2 ^( P+ i)... fa -i)g t 
X s n ... p (l-\G$ h ...fr-v)* h ..* \p + ) {l-\Gt h ... {il -i)t h ... P \(p + l) + ) ' 
/ = ji - 1 and (ji > h - 1 or j 2 = i 2 ) : , (5.37) 

/ = ji — 1 and (ji < ii — 1 and j 2 > i 2 ) : 

_ 11 t h ... h tj 1 ...j 2 |1 + ) (j2~| -^"ji...(ii-l),(i 3 +l)...j 2 ^"ii...i 2 + ) 

p=j 1 (1 I ^jl,P+l,j2,j2^P+l,W,i2,j2 |1 + ) 

(1"| ^^(p + l)...(i 1 _l),(i 2 + l)...j 2 -?t(p + l)...(i 1 _l)^ t |1 + ) 



X 



X- 



(11 ^(p+D.-.Cn-D^ib, J2 + 1) |j 2 + ) (11 ^(p+i)...^-!)^,..^ |1 + ) 

(p(p+ 1)> (11 g^(i a+ l)...jJ^ l3 - a (jl > p)] 2 ^QH.l)...(i 1 -l)g t |1 + ) 
(11 ^fe + l)... j2 ^fe J2 + 1) b + ) (11 ^fe + D.., 2 ^fe J2 + 1) l(P + 1)" 

i2 + 1 < I < 32 ■ 

(*> p> ™) + ^£u2 (*> p. n ) > 
p=j+i p=ji-i 
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and 



ji - K P < (h ~ 1) : 

(l-|jt jWa ^ 1 ,.. <a |l+) 

{(j ~ jj U1 ^(p+l)...(n-l),(»2+l)...(^-l)^(p+l)-(n-l)^ f 
X ^JW-fr-i^SfoO K Z -. 1 .) + > (11 ^(p+D-Cn-D^i^O l' + > 



(p(p+l))(l- 


^fe + i).4/-i)[^2(^2;ji,p)] 2 ^( P+ i)... (n -i)^|i + } 


(1- 






^ + i)... ( i-i)^(P'OI(p + i) + > 



p = ji - 1 and ji = h: , 
P = ji — 1 and ji < ii : 

ai^,., 2 ^...ji + ) 

((/ - 1) Q g^ J - 1 ,.. (il -i),( fa+1) „, (t -i)^ 1 „, (<1 - 1) g t |i + ) 3 

( 5 - 38 ) 

(ii ^ fe+ i)... ( /-i)^..,2^ (ii ^ + i)...(i-D^igOi - 1, o ' 

(/ + 1) < P < J2 : 

((/-i)0(i1g^... P ^ 2+ D...pg1i + ) 3 

(ii Qk^)... v ki... v \(i - 1) + ) (ii Gt {i2+ i)... P ti... ? 

(p(p+ i)) g^ fa+1) ...( t - 1) [^,p)] 2 ^ (<3+1) „, p gt |i+) 

Si...p(l~|^(i 2+ l)...(i_l)^...p |p + ) (l"|^(i 2 +l)...(i-l)^...p|(p+ 1) + ) ' 

p = h ■ 

((/-i)/)(i-|g^„, 2 ^ fe+1) .., 2 gt|i +) 3 

(11 Qt^ ) ... h ti...n l(* " 1)1 (11 ^ 2+ i)..,2^..,2 H + > 

( j2 -| g |l+> g# (ia+1) .„ (t - 1) [.T(Z, J 2 )] 2 ^ fe+1) ,., 2 g t |1+) 

S/...j2 (l"l \J2 + ) (1~| ^(i2+i)...(i-i)^i...j 2 ^ |1 + ) 

The T 3 term is given by, 

ji n—l 

= E/S2M+ E ('>")> (5-39) 

1=2 1=32 + 1 
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which depends on 



2 - /<J1 " 1: WHY) ■ 

Jl ' (0'i-i)-|^|i + ) 

/ = ji and j'2 = n '■ , 
/ = ji and j 2 <n: 

(11 ^,.. 02+ i)^ 02+ i)...n |1 + ) (11 ^ |1 + ) (1 (J2 + 1)) 



(5.40) 



(11^|(j 2 + l) + ) 



j 2 + 1 < Z < n 
Following on, the T 4 term reads, 



(l-|^ (m) ^ (m) „, w |l+) (!/)(! (/ + !)) 
(1(1 + 1)) 



n 



n-2 



/=J2 + 1 



(5.41) 



«=3 
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h + Kl<(n-2): 



n-l 



_ ^ ((/-!)/) {l-\$ L .. p $ {p+1) ... n \l + f 

p tltl ( 1 ~\^( P+ l)Jl... P \(l- + ) (l-\t {p+ l)...n$L.. P \l + ) 
(p(p+l))(l'|^,.. n [^(/, P )] 2 ^ (p+1) ... TO |l + ) 

X S L .. p {\~\tl...ntl... v \p + ) (11 $I...U$1... P \{P + + ) ' 

1 = 32 + 1- 

(l-\Q ] \(J2 + 1) + ) (1-| ^ 2+ l... P ^( P+ l)...„|l + ) 



n-l 



- v — 

p 4T +2 (1-| t (p+ l)...nt n+ l... P Q + ) (I"! ^(p + l)...n^ 2+ l...p I0'2 + I)" 



(P(P+ 1)> (11 ^ 2+ l...nTO 2 + 1,P)] 2 ^( P+ 1)...„ |1 + ) 



(1 I ^ 2 +l...n^j 2 +l...p b + ) (1 I \{P + 1) + ) ' 

I = ji and (j 2 > n - 1 or j ± = 2) : , 

/ = ji and (j 2 < n - 1 and ji > 2) : (5-42) 
g (11^,., 2 ^,., 2 |l + > 



P=J2+1 

X 



(1 I ^ji...j2-^(j2+i)...p^Ji,n,«2,P |1 + ) 



(11 ^(p+D^n^O'l - 1.P+ 1) lO'l - + ) 
J(^-l)1g|l + )(11^,.. p ^ +1) ,,,Jl + ) 3 

^( P +l)...n#(j 2 +l)...pG |1 + ) 
x (P(P+1)) 

(i1^i... (ji -i)C/S0'i-i,p) b + > 

(11 ^l-Q- 1 -l)[n&- 2 0'2 + l,p)] 2 ^( P+ l)...n |1 + ) 

<i-| ^i...o- 1 -i)^iSc/i - 1) i(p+ 1) + > 

ji n-l 

3 < I < ji : ^ 7S lh {l,P,n) + ^ 7lj lj2 (l,P,n) , 



p=l+l P=J2+1 



45 



and 



j2 + l<p<n — 1: 

(1-|^,.. J2 ^,., 2 |1 + ) 



((1-1)1) 



..p#"(p+l). 



X 



X 



X 



(l-l^... p ^ (p+1) ... n ll + ) 3 



(i-\$ {p+lh .. n u;i? 2 (i-i,p + i)\i + ) 

(p(p + l)) 

(i1^i...a-D^(/-i,p + i)b + ) 

^l...a-l)[n i ;^ 2 (^ , l - + l,p)] 2 ^( P+ l)...n |i + ) 



(l-||t M ,-i)^jS(Z-l,P+l)|(p+l) + > 
p = ji and j 2 = n: , 

p = ji and j2 < n : 

ai^,.,A.^ii + ) 

<(* -i)0 



TlJ lj2 (l,P,n) = < 



x 



x 



X 



X- 



(11^ 02+ l)...n^ 1 1 |(/-l,j2 + l)|(/-l) + ) 
(l"l^.., 2 ^Q 2+ l)...n|l + ) 3 

(r\$ (n+iy .. n u]^(i-i, ]2 + i)\i + ) 

(1-|^I0'2 + 1) + ) 



(i1^i...a-i)^(/-i,j2 + i)|(j 2 + i) + )' 

Z + 1 < p < ji - 1 : 

((/-1)/)(1-I^„, P ^ 1 „, P I1+) 3 
(l-\t,... p tl... p \(l-l) + )(l-\t,... p tl... p \l + ) 

(p(p+i))(i^$ 1 ...^ l) [Hi,p)ft l ... P \^) 



X- 



(5.43) 



ti... { i-i)ti... P \P + ) (1"| # M ,-i)#,.., \(P + 1) + ) ' 
P = ji - 1 : 

((/-i)0(i1^.., 1 -i^i.., 1 -i|i + ) 3 

S ,.. jW A' |.../_ I ^ /..../: I |(* " 1) + ) (11 |/ + ) 

(Oi - 1)1 |i + ) (il h - 1)] 2 ^!..,,-! |i + ) 



X- 



(i-||ti...(,_i)|t,.. J - 1 _i |(ii - 1) + ) Jki...{i-i)Jfci...j l -iG |i + ) 

Above we have used the following abbreviations, 

^,..* 2 |1 + ><11£ 



u}\f 2 (p,i) - ^(p+i)...(n-i),(i 2 +i)...(i-i) - ^py^^ V^) ' (5 ' ! ! 1 



as well as 



7?t 



^,., 2 |l + )(11g 
(l-|^,.. j2 ^,., 2 |l + ) 



$a...& 5 
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=1112 



(a, 6; c, d) = ^ J U2 (a, 6) + J**^ (c, d) + ^..6^ , 

6-1 6 

^(o, &) = E E ^ • ( 5 - 45 ) 



i=a m=i+l 

Next we have, 

{2 
m = 2 or m = n : - , 
2 (12)<2 3)((ra- l)n)(nl) ( 5 - 46 ) 
otherwise : - - ,„,,/,:,,, /A,. . 

9 (2n) 2 (13)((n- 1) 1) 

The unwinding of eq. (A13) gives the following contribution, 

rpi^f mm(n ^ 1 2) /C(ji + 1, j 2 ,ii,i 2 ,m)((i 1 - ljiiKjg (jg + 1)) 

^ (^nj--^ 2^ 2^ <i<i>(<2i> 

x {W&nirJiih + 1) + W^Cr, + 1, j 2 + 1) - W^r, h + 1J 2 ) - H' ^ Jr.. ,...„)). 

(5.47) 

The terms entering this expression are given by, 

(r, a, b) = (r, a, b) ± fljj ^ M) , ( 5 - 48 ) 



where 'H±*( il _ 1 jj 2 (r, ji + l,y) = ^t^-i)^^ y)- We have introduced the following func- 
tions, 



6 (a \Qh+l,h,i 2 ,32 I 1 ) I^l + Ml,i2j2 I 1 ) 



o-ni 2 /o _ n _ 1 ^(0,0,6) M 1 S JlJ2 (a,a,&) 

d (a I ^i+i,w,i2,j2 |1 ) ^ ( a I t/ji+i,n,j 2 j2 I 1 ) 

with the constraints H£j 2 (r, i 2 , x) = 0, H l ^ 2 (r,i 2 ,x) = 0, H*^(r, z, i 2 ) = and 
^|(2, a;, i 2 ) = as well as Hjlf 2 (r,x,y) = = H l -^ 2 (r,x,y) when ^Cji+i)...(ti-i),(ia+i)...j a = 
fc-j. We also need 

= -fc%^a^f%l j2 {j) {l-\g n+1M2 + ) (5.50) 

X ^Ji+l,*i,»2ja + ) + s (j'i+l)-(<i-l).(*2+l)-J'a( 1 J')) > 

X I ^'l+l,H,«2j2 + ) J 
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and 



otherwise : (1 b) (1 | G jl+ i, h ,i 2 , h \a + ) / (ab) , 



■ ■ i , • • , -i tj 1(^2 + 1)- 

j! = % x - 1 and j 2 = i2 + 1 : 



((H-l)l4 + i,n,wJ 1+ >' 



■ OH' (l-|^-|(»i-l)-> f5 53) 

ji = i 1 -2 andj 2 = ^2 : ^77^ fr i Tvk ' 1 j 



otherwise : 



S0 1 +l)...(i 1 -l),(t 2 + l)...3 2 ( 1 


^Ul+l)-^ 




1 + ) 


((n-i)- 


J1 + 1.»1.»2.J2 




^1 + 1><1><2'J2 





Finally, the large- ,2 contribution, eq. (A17), results in, 

T^ l2 (n) = (5.54) 



V V ((»1-1)«1>('2(»2 + 1)> / (ljl) 3 (lj2) 3 [j 2 Jl] 2 

-Kft&UiJ*) + ^+5u 2 0'i + 1, j 2 - 1) + Kh^fah - 1) - iC-^O'i + , 

where we set R±*j lj2 (a, b) = if = or if X a ...6 = K h ^ A2 ; otherwise R±*j lj2 (a, b) is given 
by 

™ ,„ M _ <lji> 2 <lj 2 ) 2 f {l-\$ a ..,h\^ + ) 2 ■ 

^SiyijaV"' -' O/o. o.\2/i-| rf- _ af- li+\ V / 1 I /X ^ M+\ 



(1- 


$a...b$j 2 


1 + ) 2 


(1- 


$31 $-i\...i2 


+ ) 



+ <i-|^ a ... b 4|i + >±<i-| J r a ... b 4|i + >^) 



(5.55) 



VI. CONCLUSIONS 



At the LHC, events involving large numbers of jets will play a central role in investigations 
and measurements of new physics. On the theoretical side, a proper understanding of such 
events will require NLO calculations. These calculations in turn require one-loop amplitudes 
with large numbers of hard colored final-state particles. 

In this paper we have provided new and non-trivial examples of one-loop QCD ampli- 
tudes with an arbitrary number of external gluons. We have computed these amplitudes 
using the on-shell bootstrap method [33, 34]. Previously, n-gluon amplitudes with special 
helicity configurations, with two or three color- adjacent negative-helicity legs, have been 
computed [34, 49]. In the computation of three color-adjacent negative-helicity legs, sev- 
eral theoretical issues mentioned in the introduction — non-standard channels and large 
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shift-parameter behavior — arose and were resolved in a general manner. The calculation 
presented here is in these respects simpler, in that our choice of shift eliminates these issues. 
Here we have extended the results to cover all one-loop corrections to the celebrated MHV 
amplitudes of Parke and Taylor [64] . These are amplitudes with the two negative helicities in 
arbitrary positions in the color ordering. These amplitudes are most conveniently expressed 
in a supersymmetric decomposition [7, 26] , for which many of the ingredients had been com- 
puted previously in refs. [6, 7, 18]. The ensemble of terms, including the rational parts of the 
scalar-loop contributions given here, provide an expression for this class of amplitudes that 
is compact when compared to expectations based on a brute-force diagrammatic calculation. 

These results will be of direct use in studies of multi-jet physics at the LHC. They also 
allow us to confirm the relatively mild increase in complexity of our methods as the number of 
external legs increases. For many years, it has been widely believed that the rapid growth in 
the complexity of gauge theory calculations is an intrinsic part of the perturbative expansion. 
Our construction illustrates in a non-trivial context that this is not so for the amplitudes 
presented here. 

In constructing the amplitudes we made use of a number of empirically-observed prop- 
erties. In the calculations described in this paper, we have confirmed the validity of these 
assumptions using the stringent requirements of proper symmetries and factorization in all 
channels under real momenta. For six gluons, we also confirmed agreement with the numer- 
ical results of Ellis, Giele and Zanderighi [4] and also with the analytic results of Xiao, Yang 
and Zhu [54]. 

Although analytic results for the rational parts of the two remaining, Af = non-MHV 
six-gluon amplitudes have been obtained recently [54], it is still of interest to use the recursive 
bootstrap to construct them from the known cut-containing parts [25]. One reason to do 
so is in order to rearrange the rational parts, so that their spurious singularities are more 
manifestly cancelled against those of the cut parts. A second reason would be to see whether 
a more compact form can be obtained, in the interest of faster numerical evaluation of the 
amplitude. 

We may contrast the speed of numerical evaluation using our expressions with that of the 
semi-numerical approach of refs. [2, 3]. The semi-numerical computation of the complete 
six-gluon amplitude takes 9 seconds to evaluate one J\f = helicity configuration on a 
2.8 GHz Pentium processor [4]. There are 64 helicity configurations in total. (Subsequent 
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helicity configurations will not take as long as 9 seconds, though.) We have implemented the 
M = 4, M = 1 and M = components of the six-gluon amplitude in C++, except for the 

J\f = components of the non-MHV helicity configurations — ( — I 1 h), ( 1 h+) 

and permutations (14 in all, out of 64). The terms we have implemented so far require 30 
milliseconds in total to evaluate to 9 significant digits on a 2 GHz Xeon processor, for all 64 
helicity configurations, which is considerably faster than the semi-numerical approach. 

Besides the issue of speed, there is also the question of numerical instability due to 
round-off error, near physical and unphysical kinematic singularities. Singularities that 
cancel within the Lj and Lsj functions can easily be patched, using Taylor expansions, 
in the numerical implementation of those functions. Experience with NLO programs for 
multi-jet production indicates that the relatively mild spurious singularities remaining in 
our amplitudes will not cause significant numerical difficulties. As described in ref. [65], 
the sizes of regions of numerical instability depend heavily on the powers to which singular 
denominator factors are raised. To assess the potential for numerical instabilities arising 
near the remaining spurious singularities, we consider the examples of e + e~ — > 4 jets and 
pp — > W, Z + 2 jets. The relevant one-loop amplitudes, in the form presented in ref. [21], 
have similarly mild spurious singularities, having been obtained with an early version of the 
on-shell bootstrap. In ref. [30], these amplitudes were implemented in an NLO program 
for e + e~ annihilation into four jets. No numerical difficulties were encountered because of 
the tiny size of the unstable regions. Similarly, no numerical difficulties have arisen [66] 
in the implementation of these amplitudes in the more general crossed kinematics arising 
for pp — > W, Z + 2 jet production [32]. These results suggest that jet programs using the 
amplitudes in this paper will also be free of significant complications arising from round-off 
error. 

There are a number of open issues that would be important to address. For example, it 
would be useful to have a first-principles derivation of the complex factorization properties, 
as well as of the behavior of loop amplitudes at large values of the shift parameter. In this 
regard, recent papers [67] linking tree-level on-shell recursion with gauge-theory Lagrangians 
in particular gauges may prove useful. The unitarity method with D-dimensional cuts [8] 
may also assist in the formal understanding of properties of gauge-theory amplitudes. To 
apply our techniques to non-MHV amplitudes, one must find shifts that avoid non-standard 
factorization channels [34]. Also, one should construct a satisfactory cut-completion, free of 
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spurious singularities in the shift parameter. This construction is more intricate than for 
the MHV amplitudes presented here. 

It would also be very important to apply the on-shell bootstrap to processes involving 
external vector bosons and quarks. (For it to be applicable to cases with massive quarks, one 
would need to first extend the methods to allow massive particles in the loop.) Such processes 
are of crucial importance for understanding backgrounds to new physics in supersymmetric 
and other extensions of the Standard Model. 

With a set of one-loop multi-parton matrix elements in hand, one can proceed to construct 
a numerical program for NLO differential cross sections. Although the construction of such 
programs is non-trivial, it is well-understood, and very general formalisms are known [68]. 

Another interesting open problem concerns the twistor-space properties [39] of the loop 
amplitudes. There have already been some studies of these properties [69, 70, 71]. In 
particular, the coefficients of box integrals have a surprisingly simple twistor-space structure, 
exhibiting delta-function support on intersecting lines, as described in some detail for the 
case of Af = 4 super- Yang-Mills theory [20, 23, 36, 72]. It would be interesting to map out 
the twistor-space properties of the complete amplitudes, especially in non-sup ersymmetric 
theories, using the all-n expressions obtained here and elsewhere [7, 18, 33, 34, 48, 60, 61]. 
Using the twistor-space structure as a guide, it may, for example, be possible to construct a 
set of loop-level MHV vertices, incorporating both rational and cut terms, in analogy with 
the tree-level construction of amplitudes with generic helicities using MHV vertices [15]. 

It might also be possible to use the one-loop MHV amplitudes to obtain insight into 
the size of NLO corrections to hadron collider processes with a large number of jets. At 
tree level, one of the early applications of the MHV (Parke- Taylor) amplitudes [64, 73, 74] 
was as the basis of schemes for estimating multi-jet rates at hadron colliders, in advance of 
the availability of exact matrix elements for all helicity configurations. One approach was 
to simply assume that the non-MHV amplitudes were the same as the MHV amplitudes, 
and multiply the MHV terms in the cross section by a simple combinatoric factor [75]. 
Subsequently, a procedure known as 'infrared reduction', designed to match the known 
collinear behavior, was applied to the case of four-jet production, and gave results quite 
similar to the exact matrix elements [76]. It would be interesting to see whether similar 
procedures could be applied sensibly to NLO computations as well, making use of the MHV 
loop amplitudes reported here, as well as similar ones containing quarks, should they become 
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available. 

We anticipate that the on-shell unitarity-bootstrap approach will have widespread appli- 
cations to computing the higher-multiplicity amplitudes required for next-to-leading order 
computations of phenomenological interest at the Large Hadron Collider. 
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APPENDIX A: COMPLETED-CUT TERMS OF M = MHV AMPLITUDES 



In this appendix we collect the cut parts for the Af = one-loop amplitudes with two 
negative-helicity gluons, which were obtained in ref. [18], and complete them in a convenient 
way. We also give the rational parts of these completed-cut terms in a form convenient for 
computing the overlap terms. Converting the expressions in ref. [18] to a notation similar 
to that used for the Af = 1 amplitudes in ref. [7] , and adding suitable rational terms via Lj 
functions, we have the completed-cut parts of the MHV amplitude, 
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In performing the cut completion, we have rearranged the cuts somewhat before introducing 
the Lj functions, which automatically remove the spurious singularities. The coefficients 
appearing in the cuts are, 
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The quantities b™ J2 and c™ a also appear in the expression for the Af = 1 supersymmetric 
amplitude [7]. 

The sums over the Lj triangle functions run over the ranges corresponding to all config- 
urations in fig. 11(b) with two massive corners, each containing a negative-helicity leg, 



{l,2,...,m-2}, j = m+l, 

Xj = < {1.2 //; lj. m + l<j<n, 

{2,..., m-l}, j = n, 



(A7) 
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FIG. 11: Kinematic configuration of the box and triangle functions appearing in the completed-cut 
contribution given in eq. (Al). 



Xj 



{m,m + l,...,n-l}, j = 2, 

{m, m + 1, . . . , n}, 2 < j < m — 1 , (A8) 
{m + 1, m + 2, . . . , n}, j = m — 1 . 



For Xj when m = n — 1 and j = n, and for \j when m = 3 and j = 2, the kinematics 
degenerates and the intersection of the first and third conditions in eqs. (A7) and (A8) 
should be used. The K functions correspond to cases where the two-external-mass triangles 
degenerate to bubbles. For adjacent negative helicities, there are no boxes, and the remaining 
double sum collapses to a single sum, leaving the simpler expression given in ref. [7]. 
The box function is defined by 



M„( Sl ,s 2 ;P 2 ,0 2 ) = Li 2 (l 



sis 2 / V si / V s 2 ' ^ Si 

-^-s-ms- (A9 » 

The box function is equal to minus the two-mass easy box function of ref. [21], 
M Q (s,t,P 2 ,Q 2 ) = -Ls^ e (s,£,P 2 ,(3 2 ), corresponding to a D = 6 box integral. An al- 
ternative expression for this box integral may be found in ref. [17]. 

The Lj functions [26] are given in eq. (4.6). In the above we have replaced L x with a 
somewhat different function, modified so as to better respect the symmetry properties of 
the amplitude, 

£i(r) = to(r ) (, l /r)/ 2 _ (A10) 
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In addition, we have the bubble function, 
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By using the Li and L 2 functions in the completed-cut expression (Al) instead of loga- 
rithms, we eliminate all of the spurious poles in Cff =0 that develop z dependence under the 
[l,m) shift. These poles arise from differences of four-momentum invariants, for example 



leg m but not leg 1, so they are each shifted, leading to a ^-dependent denominator. 

Note that there are a host of other spurious poles in eq. (Al) that do not develop depen- 
dence on z under the [l,m) shift. We do not need to worry about cancelling their spurious 
behavior; it will happen automatically using the rational terms we construct. We discussed 
this issue for the five-gluon amplitude in section IV A. In the n-gluon case it is a bit more 
intricate. The denominator factors of (lm) are cancelled by a numerator factor of (lm) 4 
in A^ cc (1 , m) , but in any case (lm) is left invariant by the [l,m) shift, so no ^-dependent 
pole is induced. Now, the denominator factors of (ji J2) in &™j 2 are not entirely cancelled by 
numerator factors (the function M does cancel some of their singular behavior); but again, 
because neither ji nor j 2 can equal 1 or m, these factors are left invariant by the [1, m) shift. 
Finally, there are denominator factors of (aj) and ((a + 1) j) in c™ a , e™ a and f!£ a . While j 
cannot be equal to 1 or m, a or (a + 1) can be. However, the relevant terms are protected 
from z-dependent spurious poles by numerator factors of (1 j), (m j), (ma) or (m (a+ 1)). 

Since the entire amplitude must be free of spurious poles, we can now conclude that the 
rational remainder we wish to compute is free of z-dependent spurious poles. 

The rational parts of the completed-cut terms are obtained by setting all logarithms, 
polylogarithms and n 2 terms to zero in eq. (Al), 
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After some rearrangement, this expression can be written as, 
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Applying the [l,m) shift to C^ =0 (l , m) in eq. (Al) and taking the z — > oo limit, we can 
extract the value of Inf C n required by the basic formula (3.23). The result is, 
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Up to a factor of — A^ rec (l, m), this expression can also be obtained from formula (5.54) for 
T^ l2 (n) by setting ii = i 2 — m. 
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APPENDIX B: RATIONAL PARTS OF SIX-GLUON N = MHV AMPLITUDES 



There are three independent six-gluon MHV amplitudes. All others can be obtained 
from these by cyclic permutations of the external legs. The analytic form of the M = 
amplitude with adjacent negative helicities, A^=°(l,2), may be found in ref. [33]. In this 
appendix we present analytic forms for the rational remainders of the other two independent 
MHV amplitudes. Although the structure is rather intricate, such forms can be useful for 
future phenomenological studies of four-jet events. Given their mild spurious singularities 
(compared to more direct evaluations of Feynman diagrams), as discussed in the conclu- 
sions, we do not anticipate any significant complications arising from round-off error when 
constructing an NLO program. 

The amplitude A%=°(1,3) is given by 

<f °(1, 3) = c r [C 6 (l, 3) + R 6 (l, 3) ] , (Bl) 

where 6*6(1,3) is given in eq. (Al), 
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and Inf Cq(1, 3) is given in eq. (A17) with n = 6 and m = 3. After simplification, the result 
for %1,3) is, 
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[35] 3 [46] 3 (r|(2 + 4)|3-) 3 [24] 



+ 



+ 



+ 



3 (2- 1 (3 + 5) |4-> 2 [1 4] [1 6] [3 4] [4 5] 3 <5" | (2 + 3) |4") (1 6) (2 3) (5 6) [3 4] 2 s 234 
(1-| (2 + 4) |3"> (1 2>[2 3]«1-| (2 + 4) |3"> + (1 2>[2 3]) 



6(16)(2 3)(2 5)(5 6)[3 4] 2 s 23 4 



(1-|(2 + 4)|3-) 2 (12)(15)[23] 



(14)^ 



2 (I" | (2 + 3) |4") (1 6) (2 5) 2 (5 6) [3 4]s 234 3(1 2) (2 3) (3 4) (4 5) (5 6) (6 1) 

(5-|(l + 4)|6-)(14)(15)[5 6]si46 
6s 16 (3- 1 (1+4) |6") (5"| (2 + 3) |4") (23)(25)(56) 

(12)(14) 2 [23] / 3(15)(23) \ 

6(1 6) (2 3) (2 5) (3 4) (5 6) [3 4] V (13)(25) J 

\3 /0 _ 



+ 



(2-|(l + 4)|6-r (3-|(2 + 5)|4-> [2 3] [4 6] 



3 (2"| (3 + 5) |4-> 2 (5-| (1 + 4) |6"> (5"| (2 + 3) |4") (2 3) 2 [1 4][1 6] 
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, (14) 2 (15) 2 [35] / 3(12)(35) \ 

6(1 2) (16) (2 5) (3 4) (3 5) (5 6) [3 4] V (13>(25) J 

■ (14)(15)[23][56] /(14) 3 (l'| (2 + 4) |3~) (1 5) \ 

6 (1- 1 (3 + 4) |2-> (2 5) (5 6)s 234 V (2 3) (1"| (2 + 3) |4") (2 5) J 
(11(3 + 5) |4-) 4 [35] 3 



+ 



+ 



6 (2- 1 (3 + 5) |4"> (6- 1 (3 + 5) |4") (1 2)(1 6)(3 5)[3 4] 2 [4 5] 2 s 345 

(l-|(3 + 5)|4-) 4 (26)[26][35] 4 

6 (1-| (3 + 4) |5"> (11 (4 + 5) |3"> (2"| (3 + 5) \A~f (6"| (3 + 5) |4~) 2 [34][4% 345 

(l-|(2 + 6)(3 + 5)|6+) 2 (16)[26][35] 

3 (1-| (3 + 4) |5"> (1-| (4 + 5) |3"> (6"| (3 + 5) |4"> (2 6)(3 6)(5 6)s 345 
f (14) (16) [3 5] 

^2(35) (6-|(3 + 5)|4->. 

(5- 1 (2 + 3) |6"> [5 6]s 2 46 ([l 4] (1 2) (1 5) + (2 5) (1 6) [4 6]) 



x 



+ 



3(21(3 + 5)|4-) (3-|(l + 4)|6-)(5-|(2 + 3)|4-)(16)(23)(56)[16](25)[14] 

(15) 2 [2 3] 2 [5 6] 

6 (1-| (3 + 4) |2-) (5- 1 (2 + 3) |4~> (2 5)(5 6)s 234 

/ 3(l-|(2 + 4) |3-)(15)(23) _ 2(1 5) (2 6) [2 3] \ 

V { ' (l-|(2 + 3)|4-)(25) (5 6) [3 4] J 

(14) 3 (15)[2 5] 

6(1"| (3 + 4) |2-) (3- 1 (1+4) |6") (1 3)(1 6) 2 (2 3) (2 5) 2 (3 4) (5 6) 

x 6(2 + 5) |3 + ) + (1"| (5 + 6)2 |3 + ))((1 3) (2 5) - 3(1 5)(2 3)) 



+ 



3(12)(15)(23)(35)[2 5] 

(14) 3 (3 6)[2 6] 



6 (1-| (3 + 4) |5"> (3- 1 (1+4) |2"> (1 2)(1 3) (2 6) 2 (3 4) (3 5) 2 (5 6) 
x (2((1-| 2(5 + 6) |3 + ) + (r| (2 + 6)5 |3 + ))(1 3) (2 6) 

- 3(2 2(5 + 6) |3+) + 65 |3+))(1 6) (2 3)) 

(1-1(4 + 5) |3-)(15) 2 [35] 

2 (I" | (3 + 5) |4") (1 2) (1 6) (2 3) (2 5) (3 5) (5 6) [3 4] S345 
x ^l-|(2 + 6)(3 + 5) l 2 + )(12)(35) _ (1 _ |(3 + 4)|5 . )(13)(2 ^ 



+ 



2-|(l + 4)|6-) S 2 
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3s 16 (2- 1 (3 + 5) |4-) 2 (3-| (1+4) |6") (2 3) (2 5) (2 6) [1 4] 



x 



+ 



/ (1 6) (2 5) (2 6) [4 6] 2 [5 6] (2~ | (1 + 4) |6-) (1 2) 2 [1 4] [2 6] \ 

V (5-|(2 + 3)|4-> (5-|(2 + 3)|6-> J 

(2-|(l + 4)|6-) 2 (12)[2 6]s 14 6 

2s 16 (2- 1 (3 + 5) |4-) (3- 1 (1 + 4) |6") (5"| (1 + 4) |6") (2 3) (2 5) (2 6) (5 6) [4 6] 
/ (5-|(l + 4)|6-)(15)(26) _ (6-|(l + 4)|6-) (16)(25) (14) (5 6) [4 6] ' 

V (2 5) (2 6) 3 
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, (l-|(2 + 6)(3 + 5)|6+) 2 (16)[26][35] 

+ (l-|(3 + 4)|5-)(l-|(3 + 5)|4-)(l-|(4 + 5)|3-)(6-|(3 + 5)|4-)(23)(26)(35)(36)(56) 
1 / <1-| (3 + 4) |5~> <13><2 6><3 5> <1~| (2 + 6)(3 + 5) |2+> <1 2><3 6> \ 



x 



2s 3 45 V (3 6) (2 6) ) 

(14) 3 (12)(23) 



(34)(56)(6 1) 2 (13)(35) 2 (3 + 4) |5") (3"| (1+4) |6") 

'1 (13) 2 (35) 2 (16) 2 [25] 2 [26] 2 1 (1 3) (1 5) 2 (3 5) 2 [2 5] 3 



x 



3 (12)(23) (1-| (3 + 4) |2-) (3~| (1 + 4) |2") 3 (2 3) (2 5) (3 + 4) \2~ 
1 (13)(3 6) 2 (16) 2 [2 6] 3 
3(12)(26) <3-| (1 + 4) |2") 

+ I <1 2) <2 <1 5) <2 6) (< 1 " I + 6)5 | 3+ > + <1~ | 6(2 + 5) |3^ 



x (jl 6) (3 6) (2 5) [2 6] + (15) (3 5) (2 6) [2 5] 

(13)(16) n 
(12)(26)/ 

(2 + 6)5|3+) (35) 2 [25] ((15)(23) + (13>(25» 



+ 
+ 



(2 5) 2 (2 3) 2 

(l~|6(2 + 5) |3+) (16) 2 [2 6] ((12)(3 6) + (13)(2 6)) 
(26) 2 (12) 2 

(12) 2 (2 3) 2 ' 

+ (1-|6(2 + 5)|3 + )(1-|62|3 + ))^ 



+ /1 ^2/^X2 (< 1- I ( 2 + l 3+ ) 25 l 3+ ) + 62 l 3+ ) 25 l 3+ ) 



(B7) 
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